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EXECUTIVE  SUMMARY 

A.  OBJECTIVE:  This  research  was  conducted  to  extend  our  knowledge  of  the  chemical  and 

physical  processes  which  control  the  eventual  fate  of  vapor-phase  chlorotiifluoroethylene  (CTFE) 
in  the  environment.  The  studies  were  conducted  in  a  controUed-environment  chamber  designed  to 
simulate  various  conditions  which  would  be  representative  of  those  which  CTFE  vapors  would 
encounter  in  operational  settings. 

B.  BACKGROUND:  Hydraulic  fluids  are  an  important  constituent  of  many  modem  military 
weapons  systems.  The  petroleum-based  fluids  first  used  in  hydraulic  systems  were  responsible 
for  a  significant  number  of  aircraft  fires  because  of  their  flammability.  Since  the  early  seventies, 
researchers  at  several  U.S.  military  installations  have  sought  to  develop  a  iww  class  of  hydraulic 
fluids  which  would  have  increased  performance  capabilities  and  be  nonflammable.  Several 
interim,  fire-resistant,  fluids  were  developed,  but  the  ultimate  goal  remained  an  entirely 
nonflammable  hydraulic  fluid. 

Researchers  at  the  Air  Force  Materials  Laboratory  at  Wright-Patterson  AFB  began  a  program 
in  1975  to  develop  a  new  hydraulic  fluid  with  the  required  properties.  Of  the  many  potential 
candidates,  chlorotrifiuoroethylene  (CTFE)  was  eveohially  selected  because  it  is  nonflammable  and 
has  the  required  physical  properties.  The  initial  positive  evaluation  of  CTFE  as  a  replacement 
hydraulic  fluid  has  been  tempered  by  subsequent  toxicological  studies  of  this  material.  Exposure 
studies  at  the  Toxic  Hazards  Division.  Harry  G.  Armstrong  Aerospace  Medical  Research 
Laboratory,  Wright-Patterson  AFB.  led  to  the  conclusion  that  CTFE  vapors  could  have  measurable 
toxic  effects  under  some  circumstances.  This  finding  provitted  the  impetus  for  RAVC  studies  to 
determine  the  expected  atmospheric  lifetime  of  CTFE  vapors  in  an  attempt  to  provide  additional 
data  for  an  overall  risk  assessment  of  this  material. 

C.  SCOPE:  The  study  began  with  a  determination  of  the  stability  of  CTFE  vapors  under  inert 
conditions  to  establish  its  inherent  interactions  with  the  controlled-environment  chamber .  Then  the 
effects  of  adding  oxygen,  moisture,  and  ozone  were  determined.  The  effects  of  adding  additional 
reactive  surfaces  to  the  chamber  were  also  evaluated.  Eventual  goals  were  to  modify  the  chamber 
to  include  a  powerful  light  source  which  would  enable  the  evaluation  of  surface-enhanced 
photochemical  degradation  studies.  However,  the  reassignment  of  the  principal  investigator  and 
his  assistant  led  to  a  conclusion  of  the  research  before  photochemical  studies  could  be  undertaken. 
This  report  summarizes  the  research  which  was  completed. 
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D.  METHODOLOGY :  The  research  was  conducted  in  the  Environics  Laboratory  controlled- 
environment  chamber.  CTFE  was  obtained  from  Wright-Patterson  AFB.  Species  concentrations 
in  the  chamber  were  monitored  principally  by  Fourier  transform  infrared  (FT-IR)  spectroscopy; 
some  measurements  were  also  taken  with  a  residual  gas  analyzer. 

E.  TEST  DESCRIPTION:  The  study  was  done  in  progressive  steps.  First  the  stability/reactivity 
of  CTFE  in  the  empty  chamber  in  an  inert  atmosphere  was  determined.  Then  the  effects  of  adding 
oxygen,  moisture,  and  ozone  were  measured.  Finally,  20  corroded  aluminum  plates  were  added  to 
the  chamber  increasing  the  interior  surface  area  by  about  a  factor  of  10.  CTFE  stability/reactivity 
was  then  reevaluated  with  these  plates  present 

F.  RESULTS:  Because  it  is  a  fully  halogenated  species,  CTFE  is  not  expected  to  undergo 
reactions  in  the  open  troposphere  and  its  reactivity  on  surfaces  was  expected  to  be  low.  The  results 
of  the  studies  completed  to  date  bear  out  these  expectations.  CTFE  did  not  react  in  the  presence  of 
oxygen  or  ozone  in  the  presence  or  absence  of  added  surfaces  according  to  the  results  of  this 
study. 

G.  CONCLUSIONS:  CTFE  appears  to  be  strongly  adsorbed  onto  the  surfaces  it  contacts.  It  is 
partially  desorbed  by  water  vapor.  It  does  not  appear  to  react  on  surfaces  or  in  the  atmosphere  in 
the  absence  of  sunlight.  The  strong  sorption  tendency  and  the  likelihood  of  slow  desorption  lead 
to  some  concern  about  operational  settings  where  small  spills  over  a  long  period  of  time  could 
produce  enough  vapor-phase  CTFE  to  cause  exposure  problems. 

H.  RECOMMENDATIONS:  Studies  from  the  literature  suggest  that  CTFE  would  likely  undergo 
photoassisted  degradation  reactions  on  the  soil  and  perhaps  other  surfaces.  Because  of  this 
likelihood,  eventual  project  goals  were  to  modify  the  chamber  to  include  a  powerful  light  source 
which  would  enable  the  evaluation  of  surface-enhanced  photochemical  degradation  studies. 
However,  the  reassignment  of  the  principal  investigator  and  his  assistant  led  to  a  conclusion  of  the 
research  before  photochemical  studies  could  be  undertaken.  If  possible,  the  chamber  modifications 
should  be  completed  and  the  photodecomposition  experiments  undertaken. 
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SECTION  I 


INTRODUCTION 


A.  OBJECTIVE 

This  research  was  conducted  to  extend  our  knowledge  of  the  chemical  and  physical  processes 
which  control  the  eventual  fate  of  vapor-phase  chlorotrifluoroethylene  (CTFE)  in  the  environment. 
The  studies  were  conducted  in  a  controlled-environment  chamber  designed  to  simulate  various 
conditions  which  would  be  representative  of  those  which  CTFE  vapors  would  encounter  in 
operational  settings. 

CTFE  spills  can  occur  during  both  normal  operations  and  in  accidents.  The  vapors  which 
escape  have  a  finite  lifetime,  which  depends  on  their  interactions  with  surfaces  and  in  the  open 
atmosphere.  Since  CTFE  is  a  fully  halogenated  species,  it  is  expected  to  be  unreactive  with  respect 
to  gas-phase  tropospheric  processes.  However,  its  interactions  with  surfaces  have  not  been 
characterized  and  were  the  principal  goal  of  this  effort 


B.  BACKGROUND 

Hydraulic  fluids  are  an  important  constituent  of  modem  military  weapons  systems.  The 
petroleum-based  fluids  Hrst  used  in  hydraulic  systems  were  responsible  for  a  significant  number  of 
aircraft  flres  because  of  the  flammability  of  these  fluids.  Since  the  early  seventies,  researchers  at 
U.S.  military  installations  have  sought  to  develop  a  new  class  of  hydraulic  fluids  which  would 
have  increased  performance  capabilities  and  be  nonflammable.  Several  interim,  fire-resistant, 
fluids  were  developed,  but  the  ultimate  goal  remained  an  entirely  nonflammable  '^vdruulic  fluid 
(Reference  1). 

Researchers  at  the  Air  Force  Materials  Laboratory  at  Wright-Patterson  AFB  began  a  program 
in  1975  to  develop  a  new  hydraulic  fluid  with  the  required  properties.  Of  the  many  potential 
candidates,  four  classes  of  fluids  were  able  to  meet  the  initial  performance  requirements: 
chlorofluorocarbons,  fluoroalkylethers,  linear  perfluoroalkylethers,  and  perfluoroalkylether 
triazines.  Of  these,  CTFE  was  eventually  selected  (Reference  2).  The  choice  has  also  been 
analyzed  and  documented  by  Mahony  (Reference  3).  A  brief  developmental  history  and  summary 
of  properties  are  provided  by  a  commercial  vendor  (Reference  4). 
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The  initial  positive  evaluation  of  CTFE  as  a  replacement  hydraulic  fluid  has  been  flagged  by 
subsequent  toxicological  studies  of  this  material.  Exposure  studies  at  the  Toxic  Hazards  Division, 
Harry  G.  ArmsU'ong  Aerospace  Medical  Research  Laboratory,  Wright-Patterson  AFB,  (References 
5-8)  led  to  the  conclusion  that  CTFE  vapors  could  have  measurable  toxic  effects  under  some 
circumstances.  This  finding  provided  the  impetus  for  RAVC  studies  to  determine  the  expected 
atmospheric  lifetime  of  CTFE  vapors  in  an  attempt  to  provide  additional  data  for  an  overall  risk 
assessment  of  this  material. 

C.  SCOPE 

The  study  began  with  a  determination  of  the  stability  of  CTFE  vapors  under  inert  conditions  to 
establish  its  inherent  interactions  with  the  controUed-environment  chamber  (References  9  and  10). 
Then  the  effects  of  adding  oxygen,  moisture,  and  ozone  were  determined.  The  effects  of  adding 
additional  reactive  surfaces  to  the  chamber  were  also  evaluated.  Eventual  goals  were  to  modify  the 
chamber  to  include  a  powerful  light  source  which  would  enable  the  evaluation  of  surface-enhanced 
photochemical  degradation  studies  (References  11  and  12).  However,  the  reassignment  of  the 
principal  investigator  led  to  a  conclusion  of  the  research  before  photochemical  studies  could  be 
undertaken.  This  report  summarizes  the  research  which  was  completed. 
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SECTION  II 


EXPERIMENTAL  PROCEDURES 
A.  SPHERICAL  REACTION  CHAMBER 

To  ensure  maximum  control  over  experimental  variables  in  this  study,  a  spherical,  stainless 
steel  reaction  chamber  was  employed.  The  design  and  construction  details  for  this  chamber  have 
been  previously  described  (References  9  and  10)  and  wOl  only  be  briefly  outlined  here. 

The  primary  design  consideration  for  this  chamber  was  to  minimize  surface  interactions.  To 
this  end,  a  spherical  design  was  adopted  giving  the  lowest  surface-to-volume  ratio  of  any  simple 
geometrical  solid.  The  chamber  was  also  coated  with  a  10-mil  thick  layer  of  FEP  Teflon®.  The 
size  of  the  chamber  was  chosen  to  give  a  reasonable  tradeoff  between  optical  path  length  and 
internal  volume.  The  choice  of  a  1-meter  diameter  sphere  gave  a  working  optical  path  length  of 
106  meters  while  maintaining  a  manageable  volume  of  about  525  liters. 

The  chamber  was  constructed  of  3/16-in  stainless  steel  hemispheres  welded  together.  Three 
12-inch  circular  ports  are  welded  to  corresponding  openings  in  the  chamber  wall.  Two  are 
diametrically  opposed  across  the  horizontal  diameter.  These  contain  mounts  for  the  long-path 
infrared  optics.  The  third  port  is  placed  at  right  angles  to  the  line  joining  the  other  two.  It  is  used 
for  introducing  surface  samples  into  the  chamber.  There  are  several  other  ports  in  the  chamber 
wall.  One  4-inch  optical  view  port  is  placed  above  each  of  the  12-inch  mirror  mounting  ports, 
which  facilitates  mirror  alignment  There  are  additional  ports  for  pressure  measuring  devices,  a 
temperature/humidity  probe,  a  stirring  fan.  and  a  high-vacuum  pumping  station. 

The  chamber  has  copper  coils  welded  to  its  outer  surface  for  temperature  control.  It  is  also 
surrounded  by  1.5  inches  of  Rubatex®  high-density  foam  insolation.  With  the  heating/cooling 
recirculation  unit  operating,  the  chamber  temperature  can  be  varied  from  6.4  *C  to  46.6  *C. 
Normal  operating  temperature  is  20  *C. 

A  control  panel  and  glass  sample  introduction  manifold  are  mounted  on  an  aluminum  frame 
which  is  attached  to  the  chamber  on  the  side  opposite  the  sample  introduction  port  The  control 
panel  provides  a  mount  for  the  temperature  gauge  and  for  two  different  pressure  gauges. 
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The  optical  system  for  the  chamber  was  designed  to  provide  the  maximum  possible 
pathlength  within  the  1 -meter  mirror  separation  constraint.  This  was  accomplished  by  using  a 
modification  of  the  basic  White  (Reference  13)  three-mirror  system  as  suggested  by  Olson 
(Reference  14).  The  chamber  optical  system  is  interfaced  to  a  Nicolet  Model  160-SX  Fourier 
transform  infrared  (FT-IR)  spectrometer.  The  optical  interface  diagram  and  other  details  of  the 
optical  system  have  been  provided  in  earlier  woilc  (References  9  and  10). 

An  additional  capability  was  added  to  the  chamber  for  this  research  effort,  a  residual  gas 
analyzer  (RGA)  with  an  atmospheric  pressure  sampling  accessory.  The  unit  is  a  Hiden  Analytical 
HAL  2/200  RGA  with  atmospheric  sampling  capability.  It  was  attached  to  the  spherical  chamber 
via  a  new  extension  flange  assembly  which  was  fabricated  by  Nor-Cal  Products  (Yreka,  CA) 
Two  connections  were  made  between  the  chamber  flange  and  the  RGA:  a  high-vacuum 
connection  through  a  2-inch  diameter  by  36-inch  long  flexible  stainless  steel  tube  with  a  2  3/4- 
inch,  copper  gasket,  high-vacuum  flanges  on  each  end,  and  an  atmospheric  pressure  sampling 
connection  through  a  2-foot  section  of  1/16-inch  stainless  steel  tubing.  The  chamber  with  the 
RGA  present  are  shown  schematically  in  Figure  1. 
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Figure  1.  Schematic  Diagram  of  Experimental  Apparatus 
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The  new  extension  flange  assembly  was  bolted  onto  the  sample  port  of  the  original  chamber 
(see  References  9  and  10  for  more  detailed  drawings  of  the  chamber)  sample  port  The  features  of 
the  new  extension  are  shown  in  Figure  2. 
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Figure  2.  Schematic  Diagram  of  the  Chamber  Sample  Port  Extension  Assembly. 


B.  MATERIALS 

1.  CTEE. 

The  CTFE  used  in  this  study  was  obtained  from  the  Toxic  Hazards  Division,  Armstrong 
Aerospace  Medical  Research  Laboratory,  Wright-Patterson  AFB,  OH.  It  was  transferred  to  a 
Teflon®-capped  glass  vial  and  used  without  further  puriflcation.  CTFE  is  also  known  as  3.1  oil 
because  of  its  viscosity  index.  It  is  a  mixture  of  oligomers  of  the  general  formula  Cl-(CF2CFCl)n- 
Cl.  It  has  been  shown  by  GC/MS  studies  that  CTFE  (3.1  oil  formulation)  is  a  mixture  of  two 
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primary  groups,  a  trimer  with  a  6-carbon  chain,  and  a  teU'amer  with  an  8-carbon  chain  (Reference 
7).  Based  on  these  findings,  an  approximate  molecular  weight  of  450  grams  per  mole  was  used  in 
concentration  calculations  along  with  a  density  of  1.93  grams  per  milliliter. 

2.  Methane. 

The  methane  used  in  this  study  was  Fisher  high-purity  grade  (99.99  percent). 

3.  Helium. 

The  helium  used  was  high-purity  grade  (99.995  percent). 

4.  Oxygen. 

The  oxygen  used  was  ultra-high-purity  (99.99  percent). 

5.  Water. 

The  water  used  was  laboratory  distilled  water. 

6.  Corroded  Aluminum  Plates. 

These  plates  were  prepared  and  stored  as  (tescribed  earlier  (Reference  15). 


C.  METHODS  AND  PROCEDURES 
1.  Sample  Preparation  and  Introduction. 

Samples  of  CTFE  vapor  were  prepared  by  injecting  liquid  Cl'Ft  into  the  chamber  manifold 
through  a  silicon  septum.  The  manifold  heating  tapes  were  turned  up  to  about  150  *C  to  facilitate 
evaporation  of  the  CTFE  which  was  then  flushed  into  the  chamber  with  helium.  For  runs 
conducted  under  simulated  atmospheric  conditions,  150  Torr  of  oxygen  were  added  to  610  Totr  of 
helium.  The  mixing  fan  was  operated  during  the  filling  process.  All  runs  were  made  at  760  Torr 
total  pressure  in  the  chamber  and  at  20  *C  (unless  otherwise  noted).. 
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Methane  was  added  to  some  runs  as  an  internal  standard.  It  was  added  as  a  10  cc  sample 
from  a  gas  syringe  (1  atmosphere  sample  pressure)  injected  into  the  chamber  sample  manifold 
during  the  chamber  fill  procedure. 

Water  was  added  to  the  chamber  to  create  'humid'  conditions,  i.e..  approximately  SO 
percent  relative  humidity.  This  was  accomplished  by  injecting  5  mL  of  liquid  water  into  the 
chamber  sample  manifold,  then  heating  and  evaptvating  the  water  into  the  chamber. 

Ozone  was  generated  by  using  a  commercial  production  unit  (PCI  Ozone  Corp.,  model 
HC2P-18). 

2.  FT-IR  Spectrometer  Operation. 

The  FT-IR  spectrometer  used  in  this  study  was  a  Nicolet  model  6000C  which  had  been 
upgraded  to  model  160SX  specifications.  It  includes  a  1280  processor  with  coprocessor,  256K  of 
20-bit  word  memory,  24Mbyte  hard  drive,  1  Mbyte  (8-in)  floppy  drive,  and  a  Sprite  Qnce  Read 
Many  optical  storage  accessory  (Nicolet  Maxtor  Model  RXT-800S).  Unless  otherwise  noted, 
spectra  were  run  at  1.0  cm'l  resolution  using  a  specially-fabricated  mercury-cadmium-telluride  A 
(MCT-A)  detector  (Infrared  Associates  Model  NB  FTIR).  The  number  of  coadded  scans  per 
spectrum  varied  from  16  to  256  with  128  being  the  most  common.  Spectra  were  acquired  under 
instrument  control  using  the  Nicolet  FT-IR  control  language  and  the  MACRO  commands  it 
provides. 

Molecular  species  of  interest  were  monitored  by  recording  their  absorption  spectra  as  a 
function  of  time.  A  pathlength  of  106  meters  was  used  for  most  of  the  study,  however,  this  was 
shortened  to  10  meters  during  experiments  with  high  CTFE  concentrations. 

A  background  spectrum  for  a  particular  experiment  was  obtained  by  co-adding  1024 
instrument  scans  and  storing  the  result  Subsequently  acquired  sample  spectra  were  ratioed  against 
this  stored  background  spectrum  to  give  absorbance  spectra.  A  typical  experimental  run  was 
conducted  as  follows.  Immediately  after  the  chamber  was  filled  to  begin  an  experimental  run,  the 
FT-IR  spectrometer  was  placed  under  software  control  for  automatic  data  acquisition,  processing, 
storage  and  plotting.  The  initial  spectrum  was  acquired  within  3-5  minutes  after  the  chamber  was 
filled.  Then  spectra  were  acquired  at  regular  intervals  thereafter  as  set  by  the  system  control 
software.  The  resulting  data,  including  time,  species  absorbance  values,  and  temperature,  were 
sent  to  an  open  computer  file  and  stored  in  ASCII  format  for  off  line  processing  and  plotting.  All 
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of  these  bperalions  and  calculations  were  handled  automatically  by  the  instrument's  MACRO 
software.  (See  Appendix  A  for  a  data  acquisition  flow  diagram,  and  Appendix  B-I  for  data 
acquisition  and  processing  software  program  listings.) 

3.  Ozone  Generator  Operation. 

The  ozone  generator  was  operated  by  passing  pure  oxygen  through  it  at  a  flow  rate  of  about 
SO  standard  cubic  feet  per  hour,  switching  on  the  power,  turning  up  the  'ozone  output'  control  to 
an  appropriate  setting  (usually  15-50  percent  of  full  scale)  and  allowing  the  ozone-oxygen  mixture 
to  flow  into  the  chamber.  Then  the  'ozone  output*  control  was  turned  to  zero  and  the  power  to  the 
ozone  generator  was  shut  off  while  the  oxygen  was  allowed  to  run  for  another  30-60  seconds  to 
sweep  out  any  residual  ozone  and  prevent  its  escape  into  the  lab. 

4.  RGA  Operation. 

The  RGA  system  was  operated  by  recording  a  background  spectrum  with  the  chamber 
empty,  using  the  two  most  sensitive  data  acquisition  settings  and  the  atmospheric  pressure 
sampling  inlet.  [The  high-vacuum  inlet  was  only  used  to  check  the  chamber  for  leaks,  but  was  not 
very  effective  for  reasons  which  are  not  clear  at  this  time].  Once  the  chamber  was  filled  for  an 
experimental  run,  another  sample  was  taken.  The  RGA  system  was  controlled  by  an  external 
386SX  computer  which  used  Hiden  Analytical  Ltd.  mASYST®  software  to  log  and  manipulate  the 
resulting  data. 

5.  CTFE  Infrared  Spectrum. 

The  infrared  spectnini  of  CTFE  is  dominated  by  strong  C-X  stretching  bands  (X  =  Cl,  F) 
in  the  800  -  1300  cm'l  region  (Reference  16).  Hie  spectrum,  acquired  at  1.0  cm'l,  is  shown  in 
Figure  3  with  the  principal  peak  positions  labeled.  The  most  intense  peak,  at  1205.8  cm*!  was 
used  for  quantitation  of  CTFE  during  experimental  runs. 
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Figure  3.  CTFE  Infrared  Spectrum,  1.0  cm'l  Resolution,  with  Major  Peaks  Labeled. 
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SECTION  III 


RESULTS  AND  DISCUSSION 

A.  CHARACTERIZING  CTFE  REACTIVITY  IN  THE  CHAMBER  WITH  NO 
ADDED  PLATES 

The  experiments  conducted  in  the  empty  chamber  are  discussed  in  chronological  order  in  the 
following  subsections. 

1.  CTFE  in  Diy  Helium 

Since  CTFE  is  a  fully  halogenated  species,  its  vapor-phase  reactivity  was  expected  to  be 
essentially  zero.  However,  experiments  with  other  species  in  the  chamber  have  shown  that 
adsorption  to  chamber  walls  is  a  major  factor  in  their  long-term  stability,  even  under  'inert' 
conditions.  This  also  proved  true  of  CTFE. 

Liquid  CTFE  (25-35  microliters)  was  injected  into  the  heated  chamber  manifold  and 
vaporized  then  flushed  into  the  chamber  with  dry.  pure  helium.  The  concentration  in  the  chamber 
ranged  from  5-7  parts-per-million  (ppm).  Hie  Ci'FE  concentration  was  monitored  for  12  days 
during  which  time  there  was  a  gradual  decrease  as  shown  in  Figure  4.  There  was  no  sign  of  any 
reaction  or  product  formation. 


Figure  4.  CTFE  Decay  in  Pure.  Dry  Helium  in  the  Spherical 
Chamber. 
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2.  CTFE  in  20  Percent  Oxygen  -  80  Percent  Helium 

Once  the  initial  stability  curve  had  been  qualitatively  observed,  this  run  was  conducted  to 
determine  the  potential  reactivity  of  CTFE  in  a  synthetic  air  atmosphere.  During  the  course  of  the 
experiment,  a  new  band  formed  in  the  940-1060  cm'i  region  which  appeared  to  be  a  reaction 
product  However,  its  identity  could  not  be  ascertained  immediately.  The  decay  of  CTFE  and  the 
formation  of  this  new  band  are  shown  in  Figure  5. 


Figure  5.  CTFE  in  20  Percent  Oxygen  -  80  Percent  Helium  in  the  Spherical 
Chamber.  A  Series  of  Spectra  Showing  the  Formation  of  a  New 
Species. 


After  this  initial  run  in  synthetic  air,  the  chamber  was  pumped  out  for  three  weeks.  Then  an 
experiment  was  conducted  to  see  if  there  was  any  evidence  for  CTFE  offgassing  from  the  chamber 


walls.  There  was  only  some  residual  carbon  monoxide  and  nitrous  oxide,  species  which  have 
been  observed  in  the  chamber  before  under  similar  experimental  conditions  (Reference  IS). 

Two  sequential  runs  with  CTFE  in  20  percent  oxygen  -  80  percent  helium  were  then 
conducted.  The  same  "product"  band  was  observed  in  these  runs,  strengthening  the  conclusion 
that  some  type  of  surface-catalyzed  reaction  was  occuring  in  this  system. 

Following  these  runs,  the  chamber  was  evacuated  overnight  with  the  diffusion  pump  and 
liquid  nitrogen  trap  in  line.  The  next  day  the  chambert  was  filled  with  20  percent  oxygen  and  80 
percent  helium  and  then  monitored  for  offgassing  products  for  seven  days.  The  formation  of 
CTFE  is  shown  graphically  in  Figure  6.  The  previously  observed  "product"  band  (940-1060 
cm'l)  was  observed,  along  with  small  amounts  of  CTFE,  formic  acid,  and  carbon  monoxide. 
These  results  are  shown  in  infrared  spectrum  plotted  in  Figure  7. 


Figure  6.  Chamber  Desorption  Experiment  following 

the  CTFE  in  Synthetic  Air  Runs  Showing  CTFE 
Formation. 
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Figuie  7.  Chamber  Desorptimi  Experiment  following  the 
CTFE  in  Synthetic  Air  Runs  Shows  Spectra 
Acquired  over  a  Seven-Day  Period. 


Another  set  of  experiments  was  conducted  where  CTFE  was  first  placed  in  dry  helium  and 
then  in  a  20  percent  oxygen  -  80  percent  helium  atmosphere.  In  the  absence  of  oxygen,  there  was 
only  gradual  decay  of  the  CTFE  with  no  "product”  formation.  When  oxygen  was  present,  the 
same  "product"  band  formed. 

During  the  course  of  these  experiments  in  the  empty  chamber,  the  RGA  unit  was  used  in  an 
attempt  to  provide  an  additional  analytical  monitor  for  gaseous  chamber  species.  The  sensitivity  of 
the  unit  proved  to  be  very  poor  under  the  conditions  of  a  typical  experiment  In  other  words,  with 
concentrations  in  the  ppm  range  and  pressures  near  760  Torr,  concentrations  of  minor  species  were 
too  low  to  be  seen  with  the  RGA  (future  plans  had  included  using  a  jet  separator  to  enrich  the 
samples  and  increase  the  instrument  sensitivity).  Because  of  this  problem,  one  set  of  experiments 
was  conducted  where  the  CTFE  concentration  was  increased  by  a  factor  of  10  and  methane  was 
added  as  an  internal  standard.  Typical  RGA  plots  are  shown  below  for  the  empty  chamber  (Figuie 
8)  and  for  the  chamber  with  60  ppm  of  CTFE  and  200  ppm  methane  (Figure  9). 
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Figure  7.  Chamber  Desorption  Experiment  following  the 
CTFE  in  Synthetic  Air  Runs  Shows  Spectra 
Acquired  over  a  Seven-Day  Period. 


Another  set  of  experiments  was  conducted  where  CTFE  was  first  placed  in  dry  helium  and 
then  in  a  20  percent  oxygen  -  80  percent  helium  atmosphere.  In  the  absence  of  oxygen,  there  was 
only  gradual  decay  of  the  CTFE  with  no  "product"  formation.  When  oxygen  was  present,  the 
same  "product"  band  formed. 

During  the  course  of  these  experiments  in  the  empty  chamber,  the  RGA  unit  was  used  in  an 
attempt  to  provide  an  additional  analytical  monitor  for  gaseous  chamber  species.  The  sensitivity  of 
the  unit  proved  to  be  very  poor  under  the  conditions  of  a  typical  experiment  In  other  words,  with 
concentrations  in  the  ppm  range  and  pressures  near  760  Torr,  concentrations  of  minor  species  were 
too  low  to  be  seen  with  the  RGA  (future  plans  had  included  using  a  jet  separator  to  enrich  the 
samples  and  increase  the  instrument  sensitivity).  Because  of  this  problem,  one  set  of  experiments 
was  conducted  where  the  Cl'FE  concentration  was  increased  by  a  factor  of  10  and  methane  was 
added  as  an  internal  standard.  Typical  RGA  plots  are  shown  below  for  the  empty  chamber  (Figure 
8)  and  for  the  chamber  with  60  ppm  of  CTFE  and  200  ppm  methane  (Figure  9). 


After  these  experiments,  the  chamber  was  pumped  out  for  six  weeks  and  then  filled  with  pure 
air  and  monitored  for  any  desorption  activity.  A  small  amount  of  CTFE  (about  1-2  percent  of  the 
original  amount  injected)  was  observed  to  desorb  from  the  chamber.  However,  none  of  the 
"product"  was  observed. 

3.  Ozone  Experiments 

Next,  a  set  of  experiments  was  conducted  to  determine  the  potential  for  surface-catalyzed 
reactions  between  ozone  and  CTFE.  The  first  step  in  this  procedure  was  to  see  if  ozone  was  stable 
in  the  chamber.  This  was  followed  by  additional  runs  to  see  if  there  was  any  baseline  (i.e.,  empty 
chamber)  reactivity  between  ozone  and  CTFE.  Approximately  25-50  ppm  ozone  were  introduced 
into  the  chamber  (based  on  an  absorptivity  of  9.7  cm'^atm"!,  base  e,  for  the  band  at  1055  cm'l 
[Reference  17])  in  dry  helium.  After  4  days  in  the  chamber,  the  ozone  had  decayed  to  about  one 
fourth  of  its  original  value  and  some  formic  acid  had  formed  (some  CTFE  desorption  was  also 
evident).  The  formic  acid  is  visible  in  Figure  10  which  is  the  ozone  spectrum  after  4  days  minus 
1.52  times  the  original  ozone  spectrum. 

Observation  of  the  ozone  spectra  quickly  led  to  another  interesting  conclusion:  the 
"product"  peak  in  the  earlier  CTFE  runs  in  oxygen  and  helium  was  ozone.  The  mechanism  for  its 
production  is  not  clear,  but  must  involve  a  surface-catalyzed  process  of  some  type. 


Figure  10.  Formic  Acid  from  the  Decay  of  Ozone  in  the  Spherical  Chamber 
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The  chamber  was  pumped  out  and  5  ppm  CTFE  were  added  in  750  Torr  of  dry  helium.  Then 
10  Torr  of  oxygen  were  added  with  the  PCI  ozone  generator  was  turned  on  long  enough  to 
produce  about  25  ppm  ozone.  Then  the  reaction  was  followed  for  about  8  days.  The  ozone 
decayed  in  a  nearly  linear  fashion  for  the  first  4  days,  then  abruptly  began  to  increase.  The 
increase  continued  through  the  remainder  of  the  run.  This  behavior  was  not  expected,  but  appears 
to  be  another  case  of  surface-catalyzed  ozone  production.  The  mediator  for  this  process  must  be 
CTFE,  because  in  earlier  runs  in  this  same  chamber  with  hydrazine,  this  ozone  production 
phenomona  was  never  observed. 

The  CTFE  decayed  more  rapidly  in  the  presence  of  ozone.  However,  it  showed  the  same 
general  pattern  of  decay  ~  a  faster  rate  for  the  first  few  hours,  followed  by  a  slower  rate  for  the 
remainder  of  the  run.  The  CTFE  and  ozone  time-concentration  curves  are  shown  in  Figure  1 1. 
Even  though  the  decay  rate  was  more  rapid  in  the  presence  of  ozone,  there  was  no  evidence,  in 
terms  of  discemable  product  formation,  that  any  specific  chemical  reactions  had  occurred.  The 
reason  for  the  increase  in  decay  rate  is  not  clear  at  this  time. 


TIME  (HR) 


Figure  11.  Plot  of  CTFE  and  Ozone  Concentrations  with 
Time. 


4.  High-Humidity  Experiments 


To  test  the  effects  of  humidity  on  the  decay  of  CTFE  in  the  empty  chamber,  a  run  was 
conducted  at  about  SO  percent  relative  humidity  in  helium.  Five  milliliters  of  water  were  vaporized 
into  the  chamber  and  flushed  in  with  helium.  Then  S  ppm  CTFE  were  vaporized  and  flushed  in 
with  helium.  Then  the  system  was  monitored  for  4  days.  The  CTFE  decayed  more  rapidly  than  it 
had  done  in  dry  helium;  the  decay  rate  was  about  the  same  as  when  ozone  was  present 

The  relative  decay  rates  of  the  CTFE  under  different  experimental  conditions  can  best  be 
compared  by  viewing  the  normalized  decay  curves  on  a  common  plot  This  is  shown  if  Hgure  12. 
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Figure  12.  Normalized  CTFE  Decay  Curves. 
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B.  CHARACTERIZING  CTFE  REACTIVITY  IN  THE  CHAMBER  WITH  20 
ADDED  CORRODED  ALUMINUM  PLATES 

1.  Overview  Tests 

When  the  20  corroded  aluminum  plates  (the  same  ones  which  were  used  in  the  hydrazine 
experiments  described  in  Referemce  IS)  were  added  to  the  chamber,  they  were  Hrst  isolated  and 
observed  to  be  sure  they  were  not  offgassing  any  vapors  from  previous  experiments  with 
hydrazine.  There  was  no  sign  of  any  offgassing. 

The  overview  tests  were  commenced  with  the  introduction  of  5  ppm  of  CTFE  into  the 
chamber  in  760  Ton*  of  dry  helium.  By  the  time  the  CTFE  was  introduced  (about  five  minutes),  it 
was  nearly  all  adsorbed.  The  small  amount  which  remained  decayed  away  within  three  to  four 
hours. 


To  compare  with  earlier  experiments  on  hydrazine,  the  chamber  was  pumped  down  to  500 
Torr,  then  5  milliliters  of  water  were  vaporized  and  flushed  in  with  helium  to  see  if  CTFE  would 
be  desorbed  from  the  plates  and  reappear  in  the  gas  phase.  About  25  percent  of  the  original  CTFE 
was  desorbed  and  appeared  in  the  infrared  spectrum.  Once  desorbed,  it  decayed  slowly.  This 
behavior  is  similar  to  that  displayed  by  hydrazine  and  typical  of  chemical  species  which  have  a  high 
affinity  for  surfaces. 

Finally,  the  chamber  was  pumped  to  610  Torr  and  oxygen  was  added.  There  was  no 
dramatic  increase  in  decay  rate  and  no  evidence  for  any  ozone  production. 


2.  Ozone  Experiments 

Finally,  a  run  was  conducted  where  the  chamber  was  fllled  to  400  Ton*  with  dry  helium, 
then  5  ppm  CTFE  were  vaporized  and  flushed  in  with  helium  to  740  Torr.  Then  ozone  was 
produced  as  the  chamber  was  brought  to  760  Torr  by  the  addition  of  pure  air.  The  initial  aliquot  of 
CTFE  was  rapidly  adsorbed,  so  after  about  10  minutes,  the  chamber  was  pumped  out  to  650  Torr 
and  another  5  ppm  were  injected  and  flushed  in.  The  CTFE  decayed  rapidly,  reaching  a 
concentration  of  essentially  zero  in  about  30  minutes.  The  ozone  also  decayed  rapidly  and  was 
gone  within  about  2.5  hours.  No  obvious  reaction  products  were  observed. 
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SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

A.  NO  ADDED  PLATES 

Although  the  chamber  is  Teflon®  coated,  CTFE  decay,  when  no  test  plates  were  added,  was 
still  controlled  by  surfaces  present  inside  including  the  walls,  mirrors,  and  other  external  ports  and 
openings.  This  surface  adsorption  is  essentially  irreversible  under  dry  conditions,  i.e.,  there  is 
very  little  desorption  of  CTFE  after  the  chamber  is  exposed  to  it 

When  20  percent  oxygen  is  added  to  the  helium,  the  decay  proceeds  somewhat  more  rapidly 
and  ozone  forms,  probably  from  surface-catalyzed  processes.  The  spontaneous  formation  of 
ozone  is  unusual  because  of  the  hightly  reactive  nature  of  this  species. 

When  ozone  was  specifically  generated,  it  reacted  with  CTFE  causing  the  CTFE  to  decay 
somewhat  more  rapidly  than  was  the  case  when  no  ozone  was  present  However,  no  specific 
reaction  products  could  be  detected  suggesting  that  the  ozone  serves  to  promote  surface  loss 
processes  which  remain  adsorbed. 

In  the  presence  of  large  water  concentrations  (i.e.,  runs  conducted  under  humid  conditions), 
the  decay  of  CTFE  occured  at  about  the  same  rate  as  when  ozone  was  present  under  dry 
conditions.  Again,  no  major  reaction  products  were  formed,  indicating  that  the  increased  rate  was 
restricted  to  surface  activity. 

The  observed  behavior  generally  supported  the  expectation  that  CTFE  would  be  unreactive  in 
the  vapor  phase  with  respect  to  homogeneous  processes.  Its  disappearance  in  the  chamber  can  be 
attributed  solely  to  surface  adsorption  processes,  some  of  which  are  enhanced  by  the  presence  of 
water  or  ozone. 


B.  CORRODED  ALUMINUM  PLATES 

These  plates  were  added  because  they  represented  the  most  reactive  surface  from  earlier 
experiments  with  hydrazine  (Reference  IS).  Thus  is  was  felt  that  if  any  surface-catalyzed  reactions 
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were  going  to  be  observed  with  CTFE,  they  would  probably  occur  with  these  plates.  This  proved 
not  to  be  the  case.  The  only  effect  of  the  plates  was  to  enhance  the  rate  of  the  adsorption  processes 
because  of  the  significant  increase  in  available  surface  area.  There  were  no  new  vapor-phase 
species  produced  by  surface-catalyzed  chemical  reactions.  The  only  other  difference  observed  with 
the  plates  present  was  that  no  ozone  was  produced  when  CTFE  was  introduced  into  the  chamber 
with  dry  helium  and  oxygen. 


C.  RECOMMENDATIONS  FOR  ADDITIONAL  RESEARCH 

The  experiments  summarized  in  this  report  represent  only  a  portion  of  those  which  were 
originally  planned.  Thus,  there  were  only  able  to  partially  address  the  reactivity  of  CTFE  in  the 
aunosphere.  These  studies  were  able  to  show,  however,  that  CTFE  is  probably  chemically 
unreactive  toward  surfaces  under  dark  conditions.  The  main  loss  process  appears  to  be  surface 
adsorption.  There  is  some  risk  associated  with  CTFE  desorption,  especially  after  repeated  spills  to 
the  same  area. 

The  one  area  which  requires  additional  research  is  that  of  photochemical  reactions  occurring 
on  surfaces.  There  is  strong  evidence  from  eailier  researchers  (References  1 1  and  12)  that  the 
combination  of  sunlight  and  soil  particles  (and  peritaps  other  surfaces  such  as  concrete,  asphalt, 
and  metals)  could  cause  the  photochemical  breakdown  of  CTFE. 

To  accomplish  photolysis  studies  in  the  spherical  chamber  will  require  the  installation  of  the 
Pyrex®  immersion  well  shown  schematically  in  Figure  13.  The  well  was  designed  to  accept  a  2.5 
kilowatt,  high-pressure  xenon  arc  lamp  and  provide  radiation  inside  the  chamber  in  the  290  -  700 
nanometer  range.  The  details  of  the  required  lamp  cooling  system  remain  to  be  worked  out,  but  the 
lamp  manufacturer  felt  that  the  proper  arrangements  could  be  made  (Reference  18).  There  are 
Teflon®-coated  aluminum  pans  available  for  use  in  the  chamber  which  could  serve  as  sample 
holders.  Such  tests  could  show  whether  surface-catalyzed  photolysis  processes  could  play  a  role 
in  the  environmental  fate  of  CTFE. 
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Figure  13.  Schematic  Diagram  of  Immersion  Well  for  High-Pressure  Xenon  Arc  Lamp. 
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APPENDIX  A 


FLOW  DIAGRAM  FOR  DATA 
ACQUISITION  AND  ANALYSIS 
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APPENDIX  B 


SOFTWARE  LISTING  FOR  THE  PROGRAM  USED  FOR  OVERALL 
CONTROL  OF  FT-IR  DATA  ACQUISITION  AND  STORAGE 


This  program  is  called  CT4.  It  was  written  in  the  Nicolet  MACRO  language  by 
using  the  text  editor  called  TED.  It  was  then  compiled  into  a  form  compatible  wi£  the  FT- 
IR  software  by  using  a  Nicolet  program  called  MACCRT.  The  use  of  the  text  editor  allows 
comments  to  be  included  in  a  MACRO  program  for  documentation  and  user  assistance. 

The  main  MACRO  program  (CT4)  calls  other  SUB-MACRO  programs.  There 
were  written  and  compiled  as  subsections  of  CT4.  Each  program  ends  with  the  statement 
END.  Each  new  subsection  begins  with  followed  by  the  program  name. 


!CT4 

BFN 

DFN 

VI2=0 

NSD 

OMD 

ENTER  EXCEL  FORMAT  DATE  (DETERMINE  BY  RUNNING  DATEME.FTN  IN 

NICOS) 

RTF 

OMD 

THE  CURRENT  SPECIES  TO  BE  MONITORED  ARE:  CTFE,  CH4 ,  H20,  AND 

03 

WTY=1 

SMN=1 

SIZ=1 

ITR=1 

MNT=0 

RTO=0 

OMD 

ENTER  THE  DELAY  TIME  (IN  MIN)  BETWEEN  SPECTRA 

VI3 

OMD 

PRESS  "RETURN"  TO  BEGIN 

PAU 

VI0=0 

CT7 

END 


!CT7 

FOR  AAA=1  TIL  120 
CT8 

DFN=DFN+1 

DCL=204445 

DCX 

GFN 
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TIME2  .  FTN  (FORTRAN  PROGRAM  TO  EXECUTE  THE  TIME  DELA  Y) 

FRN 

DCL=210414 

DCX 

NXT  AAA 
END 


!CT8 

SCD 

GFN 

OTEMP . FTN 
FRN 

OFN=DFN 

MOS 

RAS 

ABS 

XSP=4000 

XEP=700 

YSP=-.25 

YEP=1,25 

DSS 

XEP=0 

DFN=SFN 

VI2=VI2+1 

GFN 

EXP4 . FTN 

Frn 

GFN 

APP4 . FTN 
FRN 

DFN=OFN 

END 


(FORTRAN  PROGRAM  TO  ACQUIRE  TEMPERATURE  DATA) 


(FORTRAN  PROGRAM  TO  ACQUIRE  ABSORBANCE  DATA) 


(FORTRAN  PROGRAM  TO  STORE  ABSORBANCE  DATA) 


[NOTE:  statements  in  italics  are  comments  only,  not  a  part  of  the  actual  program] 
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APPENDIX  C 


SOFTWARE  LISTINGS  FOR  THE  FORTRAN  PROGRAM  USED  TO 
CONTROL  DATA  ACQUISITION  TIMING  FROM  1280  DATA  SYSTEM  CLOCK 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C  TIME2.ASC  PROGRAM  TO  RETRIEVE  TIME  (IN  SEC  FROM 

C  MIDNIGHT)  FROM  LOCATION  071  OF  PAGE  ZERO 

C  MEMORY  AND  USE  TO  CONTROL  THE  INTERVAL 

C  BETWEEN  SETS  OF  FT-IR  SCANS. 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C234567 


INTEGER  TIM 
INTVL=IRVAL (13765,0) 
ICOUNT=60*INTVL 
IFLAG=0 
IFLAG2=0 
ISUM=0 


C 

C  INITIALIZE  THE  TIME  COUNTER  (TIM) 
C 

A  .ZERADD  ISO  71 
A  MEMA  ISO 

TIM=IS0 


C  WRITE (2, 12)  INTVL,ICOUNT,IFLAG,TIM 

C  12  FORMAT (•  • , ’ INTVL= ‘ , 14 , 2X, ‘ ICOUNT= ’,14 , 2X, ’ IFLAG= ' 

C  ;  ,I4,2X, 'TIM=' ,I6,/) 

30  CONTINUE 

A  .ZERADD  ISO  71 
A  MEMA  ISO 

C  IF  (ISO  .GT.  TIM)  PRINT*, •IS0=’ , ISO 

IS3=IS0 

C  PRINT*, •IS3=' , IS3 

IF  (IS3  .EQ.  16)  IFLAG2=1 

C  IF  (IFLAG2  .EQ.  1)  PRINT* ,' IFLAG2= ', IFLAG2 

IF  (IFLAG2  .EQ.  1)  IS3=IS3+86400 
C  PRINT* , ■ IFLAG2= ’ , IFLAG2 , * IS3= ' , IS3 , ’ ISUM= ’ , ISUM 

IF  (IFLAG2  .EQ.  1  .AND.  (ISUM-IS3)  .LE.  0)  CALL  EXIT 
IF  (IFLAG2  .EQ.  1)  GOTO  30 

IF  (IFLAG  .EQ.  0)  ISUM=IS3+ICOUNT 

C  PRINT*, •ISUM=' , ISUM, •IS3=’ ,IS3, 
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IF  ((ISUM-IS3)  .LE.  0)  CALL  EXIT 
IFLAG=1 
GOTO  30 

END 


BATCH  FILE  USED  TO  COMPILE  THE  ABOVE  PROGRAM. 

NAME :  TIME2 . BAT 

LISTING 

RUN  FORTRAN 
TIME2.ASC:L 
RUN  RELOAD 
:B=1000-53777 
TIME2 . REL : L 

IRVAL . REL  [ ROOT , HLIB ] : L 
FORRUN.LIB  [ROOT. HLIB] :R 
NICAPL.LIB  [ROOT, HLIB] :R 
DPSYS.LIB  [ROOT, HLIB] :R 
FIOS.LIB  [ROOT, HLIB] :R 
NICSYS.LIB  [ROOT, HLIB]  :R 
/-TT;U 
/TIME2.FTN:S 
DEL  TIME2.REL 
DEL  TIME2.TMP 
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APPENDIX  D 


SOFTWARE  LISTING  FOR  THE  FORTRAN  PROGRAM  USED  TO 
READ  AND  REPORT  THE  CHAMBER  TEMPERATURE 


c 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C  PROGRAM  OTEMP.ASC 

C 

C  THIS  IS  A  PROGRAM  TO  READ  ASCII  TEMPERATURE  DATA 

C  FROM  AN  OMEGA  MODEL  680  THERMOCOUPLE  THERMOMETER 

C  INTO  THE  FT-IR  DATA  SYSTEM  AND  STORE  THE  RESULT  IN 

C  EACH  SCRATCH  FILE  STATUS  BLOCK  IN  LOCATIONS  68-71. 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

C  DEFINE  VARIABLES 

INTEGER  CHANNEL , BAUDRATE , BIT , BAUDVAL , DELAY , lOUT , I , 
IDELTA,  VALUE 
LOGICAL  VALID 

CHARACTER* 10  TMP ( 10 ) , WT , CONC 
CHARACTER*!  CHR, Tl , T2 , T3 ,T4 ,T5, T6 , T7 , T8 
C  SET  BAUDRATE  AT  1200,  CHANNEL=2  (I.E.,  CHANNEL  G) 

C  AND  CALL  SUBROUTINE  BAUDSET 
BAUDRATE=1200 
CHANNEL=2 
VALID  =.TRUE. 

C  PRINT  * , ’ CALLING  BAUDSET ' 

CALL  BAUDSET (CHANNEL,  BAUDRATE,  VALID) 

C  SEND  CHARACTER  *R"  TO  REQUEST  DATA  TRANSMISSION  FROM  THE 
C  OMEGA  UNIT 

VALUE  =82 

C  PRINT  *, 'CALLING  CHWRITE' 

CALL  CHWRITE  (CHANNEL, VALUE, VALID) 

C  RECEIVE  THE  CHARACTERS  TRANSMITTED  BY  THE  OMEGA  UNIT 
DELAY=0 
DO  10  1=1,8 

CALL  CHREAD ( CHANNEL , DELAY , CHR , VALID ) 

C  PRINT  *, ■ !=• ,I 

TMP ( I ) =CHR 

C  PRINT  *, •TEMP=' ,TMP(I) 

10  CONTINUE 

C  DO  20  1=1,8 

C  WRITE(3,100)  I,TMP(I) 

ClOO  FORMATC  '  ,  '  1=  ’  ,  I2 , 2X,  ' TEMP= '  ,  14 ) 

C20  CONTINUE 

T1=TMP(1) 
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oooooooo 


T2=TMP(2) 

T3=TMP(3) 

T4=TMP(4) 

T5=TMP(5) 

T6=TMP(6) 

T7=TMP{7) 

C  T8=TMP(8) 

WT=T1  //  T2  //  T3  //  T4  //  T5  //  T6  //  T7 
C  PRINT  *,'T=',WT,'  DEG  C* 

TEMP=REAL(WT) 

RETURN 

END 


BATCH  FILE  USED  TO  COMPILE  THE  ABOVE  PROGRAM. 

NAME:  OTEMP.BAT 

LISTING 

RUN  FORTRAN 

OTEMP . ASC : L 

RUN  FORTRAN 

BAUDSET.ASCiL 

RUN  FORTRAN 

CHWRITE.ASC:L 

RUN  FORTRAN 

CHREAD.ASC:L 

RUN  RELOAD 

:B=1000-53777 

OTEMP . REL : L 

BAUDS ET . REL : L 

CHWRITE . REL : L 

CHREAD.REL:L 

IRVAL . REL  [ ROOT , HLIB ] : L 

FORRUN.LIB  [ROOT, HLIB] : R 

NICAPL.LIB  [ROOT, HLIB] :R 

DPSYS.LIB  [ROOT, HLIB] ;R 

FIOS.LIB  [ROOT, HLIB] ;R 

NICSYS.LIB  [ROOT, HLIB] ;R 

/-TT:U 

/OTEMP. FTN:S 
DEL  OTEMP. REL 
DEL  OTEMP. TMP 
DEL  BAUDSET.REL 
DEL  BAUDSET.TMP 
DEL  CHWRITE. REL 
DEL  CHWRITE. TMP 
DEL  CHREAD.REL 
DEL  CHREAD.TMP 
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APPENDIX  E 


SOFTWARE  LISTING  FOR  THE  FORTRAN  PROGRAM  USED  TO 
PROCESS  AND  STORE  DATA  INTO  AN  ASCH  HLE  ON 
THE  NICOLET  FOR  LATER  ANALYSIS  AND  PLOTTING 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c  c 

C  EXP4.ASC  (gives  absorbance  values  for  4  species)  C 
C  THIS  PROGRAM  RETRIEVES  DATA  FROM  FT-IR  SCRATCH  FILES,  C 
C  AND  CALCULATES  BASELINE-CORRECTED  ABSORBANCE  VALUES  C 
C  FOR  ANALYTICAL  PEAKS.  THERE  IS  A  SUBROUTINE  FOR  EACH  C 
C  SPECIES  OF  INTEREST.  THERE  IS  ALSO  A  SUBROUTINES  FOR  C 
C  TEMPERATURE  DETERMINATIONS.  A  TYPICAL  LINE  OF  OUTPUT  C 
C  DATA  WOULD  LOOK  LIKE:  C 

C  C 

C  32020  08  39  15  20.5  1.034  0.0031  0.537  0.0002  C 

C  C 

C  WHERE  THE  NUMBERS  ARE  (IN  ORDER)  EXCEL  DATE,  FT-IR  TIMEC 

C  TEMPERATURE,  AND  ABSORBANCE  VALUES  FOR  CTFE,  CH4 ,  H20  C 

C  AND  03  RESPECTIVELY.  C 

C  C 

C  THE  PROGRAM  WAS  BEEN  MODIFIED  TO  PRINT  ITS  OUTPUT  TO  C 

C  A  FILE  CALLED  HOLD.DAT.  THIS  FILE  MUST  PREEXIST.'  C 

C  THE  SINGLE  LINE  OF  DATA  IN  HOLD.DAT  IS  THEN  READ  BY  C 

C  THE  PROGRAM  APP4.FTN.  APP4.FTN  THEN  READS  THE  VALUES  C 

C  FROM  THE  FILE  IROUT.DAT.  (NOTE!  IROUT.DAT  MUST  BE  C 

C  'INITIALIZED'  BY  THE  FOLLOWING  COMMAND  PRIOR  TO  THE  C 

C  START  OF  AUTOMATED  DATA  COLLECTION:  C 

C  C 

C  SX>COPY  IROUT.HLD>IROUT.DAT  C 

C  C 

C  BEFORE  APP4 . FTN  CAN  APPEND  NEW  VALUES . )  IROUT . HLD  C 

C  CONTAINS  A  HEADER  VALUE  OF  1  FOLLOWED  BY  A  LINE  OF  C 

C  ZEROES  IN  THE  PROPER  FORMAT.  APP4.FTN  THEN  REWRITES  THE  C 

C  VALUES  TO  THE  SAME  FILE  AND  AT  THE  END  WRITES  THE  VALUES  C 

C  READ  FROM  THE  FILE  HOLD.DAT.  PSUEDO  FILE  APPENDING  IS  C 

C  ACCOMPLISHED.  C 

C  C 

C  THE  INITIAL  CONTENTS  OF  HOLD.DAT  ARE  NOT  IMPORTANT  SINCE  C 

C  IT  IS  NEVER  READ  AND  IS  WRITTEN  OVER  WITH  A  SINGLE  LINE  C 

C  OF  DATA  DURING  EACH  CALL  OF  THIS  PROGRAM.  C 

C  C 

C  C 

C  PROGRAM  DIMENSIONED  FOR  1.0  CM-1  SPECTRA  C 

C  C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

REAL  Y(16384),BP1, LBP, PXA, PYA, CABS , XDIF, XDIF2 , YDIF, 

:  FAC1,FAC2,EP,SP,A1,A2,A3,A4 
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INTEGER  IDATA(512) . IDA (22528) , COUNT, DPI , DP2 , 

INTEGER  TFLAG , LDATE , NGEN, INCR 

CHARATER*2  HH,NN,SS 

CHARACTER*!  CTl , CT2 , CT3 . CT4 , C (10 ) 

CHARACTER*!  H! , H2 , N! , N2 , S! , S2 

CHARACTER* 4  CTEMP 

CHARACTER* 6 4  OUTPUT 

CHARACTER* !0  DATFIL 

INTEGER  I,  M,  NUMP,  MAXPOINTS 


SET  INITIAL  VALUES 


FLASER  =  15798.2 
NPMAX  =  22528 
SN  =  0 
NGEN  =  0 
INCR  =  1 


*********************************************************** 
READ  FTIR  PARAMETERS  FOR  DATA  RETRIEVAL: 

IDFN  =  IRVAL (14022,0) 

NSECS  =  IRVAL (14001,0) 

INODE  =  IRVAL (13004,0) 

NDPW  =  IRVAL (14000,0) 

IXSP  =  IRVAL (14025,0) 

IXEP  =  IRVAL (14026,0) 

IRFN  =  IRVAL (14030,0) 

IVI2  =  IRVAL (13764,0) 


THESE  IRVAL  FUNCTIONS  READ  THE  MEMORY  LOCATIONS 
CORRESPONDING  TO  THESE  PARAMETERS  (ALL  ARE  INTEGERS)  : 


(OCTAL)  LOCATION 
13004 

14000 

14001 
14022 

14025 

14026 


PARAMETER 
INODE  NO. 

NDP  (=INTEGER/256) 

FSZ  (=INTEGER/512) 

DFN 

XSP 

XEP 


READ  FTIR  PARAMETERS  FOR  SETTING  SUBROUTINE  FLAGS: 
(1=CTFE,  2=  CH4,  3=  H20,  4=  03,  5=  ??,  6=  ??) 

IFLAGl  =  IRVAL (13710,0) 

IFLAG2  =  IRVAL (13713,0) 

IFLAG3  =  IRVAL (13714,0) 

IFLAG4  =  IRVAL (13715,0) 
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IFLAG5  =  IRVAL (13731,0) 

IFLAG6  =  IRVAL (13750,0) 

C 

C 

C  THESE  IRVAL  FUNCTIONS  READ  THE  MEMORY  LOCATIONS 
C  CORRESPONDING 


C  TO 
C 

THESE  PARAMETERS  (ALL 

ARE  INTEGERS) 

C 

(OCTAL)  LOCATION 

PARAMETER 

C 

13710 

WTY 

C 

13713 

SMN 

C 

13714 

SIZ 

C 

13715 

ITR 

C 

13731 

MNT 

c 

13750 

RTO 

C 

Q*  1i  *  1i  h  *  -k  *  *  *  ■»  *******  *  *  h  *  *  *  *  it  *  *  *  *  -kit  ***************************  * 

C 

C  THE  FIRST  88  SECTORS  OF  THE  SCRATCH  AREA  (I.E.,  SECTORS  0- 
C  87)  CONTAIN  MISCELLANEOUS  INFORMATION.  SCRATCH  FILE  0 
C  BEGINS  AT  SECTOR  88  AND  OCCUPIES  FSZ/512  SECTORS.  ALL 
C  SUBSEQUENT  FILES  ARE  THE  SAME  SIZE  AS  FILE  0. 

C 

C  THE  LAST  SECTOR  OF  EACH  FILE  IS  THE  FILE  STATUS  BLOCK. THERE 
C  ARE  A  CONSIDERABLE  NUMBER  OF  SECTORS  BETWEEN  THE  END  OF  THE 
C  DATA  AND  THE  BEGINNING  OF  THE  FILE  STATUS  BLOCK.  IT  IS  NOT 
C  CLEAR  WHAT  THESE  EXTRA  SECTORS  ARE  USED  FOR. 

C 

Q*  **********************************************************  * 

c 

C  LOCATE  SECTORS  TO  BE  READ  USING  THE  FOLLOWING  INDEX. 

C 

IDFSEC= (IDFN+1) *NSECS+87 
C  PRINT  998,  IDFSEC 

C  998  FORMAT (•  ' , ' IDFSEC=  ',16,//) 

C 

C 

************************************************************* 

C 

C  READ  HEADER  BLOCK  INFORMATION  USING  IRTISK.  STORE  RESULTS 
C  IN  AN  ARRAY  CALLED  IDATA. 

C 

CALL  IRTISK ( IDATA, 512 , IDFSEC , INODE) 

C  DO  10  1=1,17 

C  PRINT  1000,  I,IDATA(I) 

CIOOO  FORMATC  IDATA*  ,  I2 , 2X,  •  =  *,16) 

C  10  CONTINUE 
C 

C  NOW,  FILE  STATUS  INFO  CAN  BE  DETERMINED  FROM  THE  ELEMEMTS 
C  OF  IDATA: 

C  lABS  =  ABSORBANCE  FLAG,  0  IF  ABSORBANCE 

C  IPTS  =  NO.  OF  DATA  POINTS 

C  FEXP  =  EXPONENT  FOR  FILE  SCALING 

C  FXAX  =  DATA  POINT  SPACING  ON  X  AXIS 
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c 

lABS  =  IDATA(IO) 

IF  (lABS.NE.O)  GOTO  900 
IPTS  =  IDATAdS)  *128 
FEXP  =  2** {19-IDATA(6) ) 

FXAX  =  FLASER/ (IPTS*IDATA(17) ) 

C 

C - 

c 

C  NOW  DETERMINE  THE  STARTING  AND  ENDING  POINTS  FOR  THE  FILE 
C  TO  BE  READ.  THIS  IS  DONE  BY  TWO  SEQUENTIAL  IF  -  ENDIF 
C  LOOPS. THE  FIRST  SETS  ISTP  EQUAL  TO  THE  LESSER  OF  XSP  AND 
C  XEP  (WHICH  ARE  READ  FROM  THE  FT-IR  SCRATCH  FILE  HEADER 
C  BLOCK) .  THE  SECOND  LOOP  CONVERTS  THE  XSP  AND  XEP  VALUES 
C  FROM  CM-1  TO  SCRATCH  FILE  DATA  POINT  VALUES.  THE  0.5 
C  FACTOR  WHICH  IS  ADDED  ENSURES  THAT  INTEGER  ISTP  AND  lENP 
C  VALUES  ARE  ROUNDED  CORRECTLY. 

C 

IF  (IXSP.LE.IXEP)  THEN 
ISTP  =  IXSP 
lENP  =  IXEP 
ELSE 

ISTP  =  IXEP 
lENP  =  IXSP 
ENDIF 

IF  ( IXSP. GT. IXEP)  THEN 
ISTP  =  IXEP/FXAX+0.5 
lENP  =  IXSP/FXAX+0.5 
ELSE 

ISTP  =  IXSP/FXAX+0.5 
lENP  =  IXEP/FXAX+0.5 
ENDIF 

NP  =  IENP-ISTP+1 

PRINT  850,  ISTP,IENP,NP 

850  FORMAT (■  ISTP=  * , I6, 2X, ' IENP=  ■,I6,2X, •NP=  ’ 

&,I6,/) 

CHECK  FOR  A  REASONABLE  NO.  OF  PTS.  IN  NP 

IF  (NP.GT.NPMAX)  THEN 
WRITE (2, 901)  NPMAX,NP 
ELSEIF  (NP.LT.5)  THEN 
WRITE (2, 902)  NP 
ENDIF 
C 

C  NOW  READ  DATA  FROM  SCRATCH  FILE  DFN  INTO  ARRAY  IDA 
C 

C  IN  IS  A  COUNTER.  NDONE  TAKES  THE  DATA  POINT  VALUE  OF  ISTP 
C  AND  CONVERTS  IT  TO  SECTORS.  NSKIP  TAKES  THE  STARTING  DATA 
C  POINT  VALUE  (ISTP)  AND  SUBTRACTS  AN  INTEGER  NO.  OF  DATA 
C  POINTS,  BY  SECTOR,  FROM  IT,  LEAVING  THE  POINT  IN  THE 
C  SELECTED  SECTOR  WHERE  DATA  REPRESENTING  THE  INTERVAL  ISTP  - 
C  lENP  BEGIN. 

C 
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IN  =  1 

NDONE  =  ISTP/512 
NSKIP  =  ISTP-NDONE*512 
PRINT  400,  NDONE. NSKIP 

400  FORMAT (■  NDONE=  ‘ , I8 , 2X, • NSKIP=  *,18,/) 

IF  THE  SCRATCH  FT-IR  FILE  HAD  XSP=0  OR  XEP=0,  THEN  DATA  CAN 
BE  READ  DIRECTLY,  ONCE  THE  PROPER  BEGINNING  SECTOR  IS 
SPECIFIED.  IF  XSP  OR  XEP  WERE  NOT  ZERO,  THEN  NDONE  SECTORS 
PLUS  NSKIP  POINTS  MUST  BE  SKIPPED  OVER  TO  BEGIN  READING 
DATA. 


RESET  FLAG  FOR  BASELINE  CORRECTION  ROUTINE 
IFLAG=0 

IF  (NSKIP. GT.O)  THEN 

SET  FLAG  FOR  BASELINE  CORRECTION  ROUTINE 
IFLAG=1 

CALL  IRTISKdDA,  512 , 88+IDFN*NSECS+NDONE,  INODE) 

DO  20  I=NSKIP+1,512 
IDA(IN)  =  IDA(I) 

IN  =  IN  +  1 
20  CONTINUE 

NDONE  =  NDONE  +  1 
ENDIF 

NPREM  IS  A  COUNTER  WHICH  TELLS  HOW  MANY  DATA  POINTS  TO  READ 
FROM  THE  FT-IR  SCRATCH  FILE  TO  COVER  THE  INTERVAL  SPECIFIED 
BY  ISTP  AND  lENP. 

NPREM  =  NP-IN+1 
PRINT  402,  IN, NPREM 

402  FORMAT {'  ’ , * IN=  ' , I6, 2X, •NPREM=  ’,110,/) 

IF  (NPREM. GT.O)  THEN 
CALL 

IRTISK ( IDA ( IN) , NPREM, 88+IDFN*NSECS+NDONE, INODE) 

ENDIF 

C 

C  NOW  THE  POINTS  IN  THE  INTEGER  ARRAY  IDA  ARE  CONVERTED  TO 
C  REAL  VALUES  BY  DIVIDING  EACH  ONE  BY  FEXP  (SEE  FT-IR 
C  SOFTWARE  MANUAL,  P.  18-5  FOR  AN  EXPLANATION).  NOTE  THAT 
C  ANY  TOTALLY  ABSORBING  DATA  POINTS  HAVE  THE  VALUE  -524288 
C  AND  THESE  ARE  SET  EQUAL  TO  ZERO  IN  THE  REAL  DATA  ARRAY. 

C 

DO  50  I  =  1,NP 
IDAS  =  IDA (I) 

IF  ( IDAS. EQ. -524288)  IDAS  =  0 
Y(I)  =  IDAS/FEXP 
C  PRINT  403,I,Y(I) 

C403  FORMAT (•  ', 'X (', I3 ,*)=•, F6 . 3 ) 

50  CONTINUE 
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C  DO  51  1=1980,2032 

C  WN= (I-l) *FXAX 

C  PRINT  71,WN,Y(I) 

C  71  FORMATC  '  .  ‘  WN=  *  ,  F8 . 3 , 2X,  ‘  Y=  ‘  ,  F8 . 3  ) 

C  51  CONTINUE 

C 

Q  *********************************************************** 

C  NOW  THAT  THE  DATA  POINTS  ARE  AVAILABLE,  BASELINE  CORRECTION 
C  CAN  BE  MADE.  THE  CORRECTION  WILL  BE  MADE  IN  THE  SAME 
C  MANNER  IT  WOULD  BE  DONE  BY  HAND.  A  POINT  ON  ONE  SIDE  OF 
C  THE  ANALYTICAL  PEAK  WILL  BE  CHOSEN  BY  AVERAGING  THE 
C  ABSORBANCE  VALUES  OVER  A  CERTAIN  PORTION  OF  THE  X-AXIS. 

C  FOR  SPECIES  WHICH  HAVE  HIGHLY  STRUCTURED  PEAKS,  THIS  MAY  BE 
C  A  VERY  SHORT  DISTANCE  (E.G.,  METHANE  USES  A  2  CM-1  ONE  ONE 
C  SIDE  OF  THE  ANALYTICAL  PEAK  AND  A  1.0  CM-1  INTERVAL  ON  THE 
C  OTHER.  THESE  INTERVALS  ARE  BASED  ON  EXPANDED  PLOTS  AND  ARE 
C  CHOSEN  TO  GIVE  POINTS  CLOSE  TO  THOSE  ONE  WOULD  CHOOSE  IF 
C  DOING  THE  BASELINE  CORRECTION  BY  HAND) .  THE  SAME  THING  IS 
C  DONE  ON  THE  OTHER  SIDE  OF  THE  ANALYTICAL  PEAK.  THESE 
C  POINTS  WILL  BE  DEFINED  WITHIN  THE  SUBROUTINE  FOR  EACH 
C  SEPARATE  SPECIES. 

Q************************************************************ 

C 

c 

C  TO  CONFORM  WITH  LOTUS  1-2-3  DATE  FORMAT,  READ  THE  VALUE  OF 
C  RTP,  THEN  PRINT  THIS  AS  LDATE 
C 

LDATE  =  IRVAL (13751,0) 

C  TO  AVOID  EXTRA  CHARACTERS  IN  THE  TIME  PRINTOUT, 

C  PICK  OUT  THE  CHARACTERS  NEEDED  FROM  THE  IDATA  ARRAY  AND 
C  CONCATENATE  THEM  INTO  THE  NECESSARY  OUTPUT. 

C 

DO  25  1=58,65 

IF  (I.EQ.58)  THEN 

Hl=CHAR ( IDATA ( I ) ) 

ELSEIF  (I.EQ.59)  THEN 
H2=CHAR(IDATA(I) ) 

ELSEIF  (I.EQ.61)  THEN 
N1=CHAR(IDATA(I) ) 

ELSEIF  (I.EQ.62)  THEN 
N2=CHAR(IDATA(I) ) 

ELSEIF  (I.EQ.64)  THEN 
S1=CHAR{IDATA(I) ) 

ELSE  (I.EQ.65) 

S2=CHAR(IDATA(I) ) 

ENDIF 

25  CONTINUE 

CONCATENATE  THE  APPROPRIATE  IDATA  ELEMENTS  TO  FORM  TIME 
CHARACTERS  WHICH  CAN  BE  PRINTED. 

HH  =  HI  //  H2 
NN  =  N1  //  N2 
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SS  =  SI  //  S2 


C 

C  RETRIEVE  THE  TEMPERATURE  CHARACTERS  FROM  THE  IDATA  ARRAY 
C  (THESE  CHARACTERS  ARE  WRITTEN  TO  THE  IDATA  ARRAY  DURING 
C  DATA 

C  ACQUISITION  UNDER  THE  MACRO  DAS  WITH  THE  FORTRAN  PROGRAM 
C  "OTEMP.FTN" ) 

C 

CT1=CHAR ( IDATA ( 68 ) ) 

CT2  =CHAR ( IDATA (69)) 

CT3  =CHAR { IDATA (70)) 

CT4  =CHAR ( IDATA (71)) 

NOW,  CONCATENATE  THESE  CHARACTERS  TO  FORM  THE  TEMPERATURE 
CTEMP=CT1  //  CT2  //  CT3  //  CT4 


DETERMINE  WHICH  SPECIES  ARE  TO  BE  ANALYZED. 
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IF  (IFLAGl.EQ.l)  CALL 
Al=CABS 

PRINT*, 'COMPUTED  Al * 

IF  (IFLAG2.EQ.1)  CALL 
A2=CABS 

PRINT* , ' COMPUTED  A2 ' 

IF  (IFLAG3.EQ.1)  CALL 
A3=CABS 

PRINT* , • COMPUTED  A3 ' 

IF  (IFLAG4.EQ.1)  CALL 
A4=CABS 

PRINT* , • COMPUTED  A4 ' 

IF  (IFLAGS.EQ.l)  CALL 
A5=CABS 

IF  (IFLAG6.EQ.1)  CALL 
A6=CABS 


CTFE(Y,FXAX,CABS) 

CH4 (Y,FXAX,CABS) 

H20(Y,FXAX,CABS) 

03 (Y,FXAX,CABS) 

??  (Y,FXAX,CABS) 
??  {Y,FXAX,CABS) 


PRINT  155 , LDATE, HH, NN, SS , CTEMP, Al , A2 , A3 , A4 
PRINT*, 'COMPLETED  PRINT  TO  SCREEN! ! ' 

OPEN ( 11 , FILE= ' HOLD . DAT ' , STATUS= ' OLD ' ) 

WRITE (11, 150)  LDATE,  HH,NN, SS, CTEMP, Al , A2 , A3 , A4 
PRINT*, 'COMPLETED  WRITE  TO  HOLD.DAT! ! ' 

150  F0RMAT(I5,2X,3 (A2,2X) , A4, 2X, F6 . 3 , 2X, F6 . 3 , 2X, F6 . 3 , 
&  2X,F6.3) 

155  F0RMAT(1X,I5,2X,3 (A2,2X) , A4, 4 (2X, F6 . 3 ) ) 


CLOSE (11) 

GOTO  999 

FORMAT  STATEMENTS 

900  WRITE (2, 990)  IDATA (10) 

901  FORMAT ('  ',/,'  ONLY  DIMENSIONED  TO  ',I6,'  PTS ' , '  NOT  ' 
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&,I6, /.  ) 

902  FORMAT (•  * , 'WHY  ARE  THERE  ONLY  ',16, •  PTS . 
REQUESTED? ' , / , ) 

990  FORMAT {•  ' ,11,'  ABSORBANCE  FLAG  =  ■,I2,2X, ‘NOT  AN  ' 
& ,  •  ABSORBANCE  FILE ',!!') 

999  CONTINUE 
CALL  EXIT 
END 


BATCH  FILE  USED  TO  COMPILE  THE  ABOVE  PROGRAM. 

NAME:  EXP4.BAT 

LISTING 

RUN  FORTRAN 
EXP4.ASC:L 
RUN  FORTRAN 
CTFE.ASC:L 
RUN  FORTRAN 
CH4 . ASC : L 
RUN  FORTRAN 
NH3.ASC:L 
RUN  FORTRAN 
H20.ASC:L 
RUN  RELOAD 
:B=1000-53777 
EXP4 . REL ; L 
CTFE.REL:L 
CH4 . REL : L 
NH3 . REL : L 
H20.REL:L 

IRVAL . REL  [ ROOT , HLIB ] ; L 

FORRUN .LIB  [ ROOT , HLIB ] : R 

NICAPL.LIB  [ROOT, HLIB] :R 

DPSYS . LIB  [ROOT, HLIB] : R 

FIOS.LIB  [ROOT, HLIB] :R 

NICSYS.LIB  [ROOT, HLIB] :R 

/-TT:U 

/EXP4.FTN:S 

DEL  EXP4.REL 

DEL  EXP4.TMP 

DEL  CTFE.REL 

DEL  CTFE.TMP 

DEL  CH4.REL 

DEL  CH4.TMP 

DEL  NH3.REL 

DEL  NH3.TMP 

DEL  H20.REL 

DEL  H20.TMP 
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APPENDIX  F 


SOFTWARE  LISTING  FOR  THE  FORTRAN  PROGRAM  USED  TO 
STORE  DATA  IN  AN  ASCU  FILE  ON  THE  NICOLET  FOR 
LATER  ANALYSIS  AND  PLOTTING 


cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C  APP4.ASC  (stores  absorbance  values  for  4  species) 

c 

C  THIS  PROGRAM  WAS  DESIGNED  TO  READ  THE  OUTPUT  FILE 
C  FROM  THE  PROGRAM  EXP4.FTN.  THE  NAME  OF  THE  EXP4.FTN 
C  OUTPUT  FILE  IS  HOLD.DAT.  IT  CONTAINS  A  SINGLE  LINE 
C  OF  DATA  WHICH  CONSISTS  OF  DATE,  TIME,  TEMPERATURE, 

C  AND  ABSORBANCE  VALUES  FOR  4  SPECIES 
C 

C  THIS  PROGRAM  THEN  READS  THE  PREEXISTING  DATA  FILE 
C  IROUT.DAT.  THIS  FILE  CONTAINS  ALL  OF  THE  PREVIOUS 
C  DATA  WRITTEN  TO  THE  FILE  HOLD.DAT.  ONCE  ALL  OF  THE 
C  DATA  IS  READ,  THE  FILE  IS  THEN  OVERWRITTEN  WITH  THE 
C  ADDITION  OF  THE  SINGLE  LINE  OF  DATA  JUST  READ  FROM 
C  HOLD.DAT.  THE  RESULT  IS  A  ROUND  ABOUT  WAY  OF 
C  APPENDING  DATA  TO  A  FILE. 

C 

C  AT  THE  TOP  OF  THE  FILE  IROUT.DAT  IS  A  HEADER  WHICH 
C  CONSISTS  OF  THE  NUMBER  OF  ENTRIES  THUS  FAR  IN  THE 
C  FILE  IROUT.DAT.  WHEN  STARTING  A  RUN,  THIS  FILE  MUST 
C  CONTAIN  AT  LEAST  A  SINGLE  LINE  OF  DUMMY  VALUES  IN 
C  THE  PROPER  FORMAT. 

C 

C  A  MAXIMUM  OF  400  SAMPLES  CAN  BE  HANDLED  BY  THIS 
C  PROGRAM. 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

INTEGER  HH,NN,SS 
REAL  CTEMP 
INTEGER  HUMID 
REAL  Al , A2 , A3 , A4 

INTEGER  I,  M,  NUMP,  MAXPOINTS,  LDATE 
PARAMETER (MAXPOINTS=400) 

PARAMETER (M= 5 ) 

REAL  RDATl (MAXPOINTS) ,  RDAT2 (MAXPOINTS ) 

REAL  RDAT3 (MAXPOINTS) ,  RDAT4 (MAXPOINTS ) 

INTEGER  CHARDAT (MAXPOINTS,  M) 

INTEGER  LDATEVECT (MAXPOINTS ) 

REAL  CTEMPVECT (MAXPOINTS) 

GET  THE  NEW  CALCULATED  DATA  FROM  DISK 
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PRINT*, 'APP.FTN  -  BEFORE  HOLD.DAT  OPENED' 

OPEN ( 11 ,  FILE= ' HOLD.DAT  * ,  STATUS= ' OLD ' ) 

READdl,  155)  LDATE,  i!H.  NN.  SS,  CTEMP,  A1 ,  A2  ,  A3,  A4 
CLOSE (11) 

OPEN  AND  READ  FILE  OF  PAST  VALUES 

PRINT*, 'APP.FTN  —  BEFORE  IROUT.DAT  OPEN' 

OPEN (11,  FILE= ' IROUT . DAT ' , STATUS= ‘ OLD ' ) 

READ (11,  140)  NUMP 

PRINT*, 'APP.FTN  --  VALUE  FOR  NUMP=  ',NUMP 
IF  (NUMP.GT.O)  THEN 
DO  860,  1=1,  NUMP 

READdl,  155)  LDATEVECT(I)  ,  CHARDAT  ( 1 , 1 )  , 

&CHARDAT ( 1 , 2 ) ,  CHARDAT ( 1 , 3 )  ,  CTEMPVECT ( I ) ,  RDATl(I), 
&RDAT2  (I) ,  RDAT3 (I) ,  RDAT4 (I) 

860  CONTINUE 
ENDIF 
CLOSE (11) 

REOPEN  THE  SAME  FILE  AND  WRITE  DATA  +  NEW  DATA 


NUMP  =  NUMP+1 

OPEN  ( 11 , FILE*  * IROUT . DAT ' , STATUS=  *  OLD ' ) 

WRITE (11,  140)  NUMP 
IF  (NUMP.GT.l)  THEN 
DO  870,  1=1,  NUMP-1 

WRITEdl,  155)  LDATEVECT(I)  ,  CHARDAT(I,  1), 
&CHARDAT(I,2) ,  CHARDAT(I,3)  ,  CTEMPVECT ( I ), RDATl ( I ) , 
&RDAT2 ( I ) ,  RDAT3 ( I ) ,  RDAT4 ( I ) 

870  CONTINUE 
ENDIF 

WRITEdl, 155)  LDATE,  HH ,  NN ,  SS ,  CTEMP ,  Al,  A2,  A3,  A4 
140  FORMAT (15) 

155  F0RMAT(I5,2X,3 (I2,2X) , F4 . 1 , 2X, 4 (F6 . 3 , 2X) ) 

CLOSE (11) 

CALL  EXIT 
END 
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APPENDIX  G 


SOFTWARE  LISTING  FOR  THE  FORTRAN  PROGRAM  USED  TO 
CALCULATE  THE  BASELINE-CORRECTED  ABSORBANCE 
VALUE  FOR  CTFE  VAPOR  AT  1205  CM-1 


SUBROUTINE  CTFE (Y, FXAX, CABS) 

REAL  Y(22528) , BPl, BP2 , PXA, PYA,CABS , XDIF, XDIF2 . YDIF, 
1FAC1,FAC2,EP,SP, IC 
INTEGER  COUNT 
C 

C  TAKES  THE  SET  BASELINE  POINTS  AND  CALCULATES  A  BASELINE- 
C  CORRECTED  ABSORBANCE  VALUE  FOR  THE  ANALYTICAL  PEAK  AT 
C  1205  CM-1.  THEN  PRINTS  OUT  THE  VALUE. 

C 

C 

C  DEFINE  THE  BASELINE  POINTS  FOR  THE  CTFE  ANALYTICAL  PEAK 
C 

IPT1=1309/FXAX 
IPT2=1311/FXAX 
IPT3=1080/FXAX 
IPT4bl085/FXAX 

PRINT  35,  IPT1,IPT2,IPT3,IPT4 

35  FORMAT (•  ' , ' IPTl= ' , I4 , 2X, ' IPT2= ' , 14 , 2X, ' IPT3= ’ , 14 , 2X, 
l'IPT4=' ,14,/)  C 

AVERAGE  THE  POINTS  FOR  THE  1ST  BASELINE  VALUE 

SP  =  0 
BP1=  0 
COUNT  =  0 
IFIRST=IPT1 
ILAST  =IPT2 

NOW  CALCULATE  FIRST  BASELINE  POINT 
PRINT  36,  IFIRST, ILAST 

36  FORMATC  ■  ,  '  IFIRST=  ’  ,  I4, 2X,  •  ILAST=  ’  ,  14 ,  /  ) 

DO  70  1= I FIRST, ILAST 

SP=SP+Y(I) 

COUNT=COUNT+l 

CALCULATE  WAVENUMBER  VALUES  FOR  PRINTOUT 

WN= (I-l) *FXAX 
PRINT  37,WN, Y(I) ,SP,COUNT 

37  FORMATC  •,'WN=',  F8 . 3 ,  2X,  ■Y(I)  =  *,  F8 . 3 ,  2X, 
&'SP=’ ,F8.3, 

12X, ■COUNT=' ,16) 

70  CONTINUE 


42 


BP1=SP/C0UNT 

C 

C  AVERAGE  THE  POINTS  FOR  THE  2ND  BASELINE  VALUE 
C 

EP  =  0 
BP2  =  0 
COUNT  =  0 
ISTART=IPT3 
lEND  =IPT4 
C 

C  NOW  CALCULATE  THE  SECOND  BASELINE  POINT 
C 

C  PRINT  39,  ISTART.IEND 

C  39  FORMATC  ‘  '  ISTART=  ‘  ,  14 , 2X,  *  lEND  =  ‘  ,  14 ,  /  ) 

DO  80  I=ISTART,IEND 
EP=EP+Y{I) 

COUNT=COUNT+l 

C 

C  CALCULATE  WAVENUMBER  VALUES  FOR  PRINTOUT 
C 

WN= (I-l) *FXAX 

C  PRINT  38. WN,Y(I) ,SP, COUNT 

C  38  FORMATC  '  ,  *  WN=  *  ,  F8 . 3 , 2X,  •  Y  (I)  =  •  ,  F8 . 3  , 

C  2X, 'SP=' ,F8.3, 

C  12X, •COUNT=' ,16) 

80  CONTINUE 

BP2=EP/COUNT 

C 

C  PRINT  OUT  THE  BASELINE  DETERMINING  POINTS 
C 

C  PRINT  85,  BP1,BP2 

C  85  FORMATC  ■,'BPl=  '  ,  F9 . 4 , 2X,  ’  BP2=  •,F9.4,/) 

C 

C  DETERMINE  THE  X,Y  VALUES  (CM-1  AND  ABS.)  FOR  THE 
C  ANALYTICAL  PEAK  (1205  CM-1)  BY  FINDING  THE  MAXIMUM  Y  VALUE 
C  IN  THE  INTERVAL  FOR  THIS  PEAK  WHICH  IS;  1207  CM-1  TO  1203 
C  CM-1  (OR  IN  DATA  POINT  VALUES;  2503  TO  2495) . 

C 

PXA=0 

PYA=0 

DO  82  1=2495,2503 

IF  (Y(I) .GT.Y(I-l) )  THEN 
PYA=Y(I) 

IC=I-1 

PXA=IC*FXAX 

ENDIF 

C  PRINT  299, I,Y(I) ,PYA, IC.PXA 

C299  FORMATC  '  ,  '  1=  '  ,  I4 , 2X,  *  Y=  *  ,  F8 . 3 , 2X,  '  PY/.^  •  ,  F8 . 3 , 2X, 

C  l‘IC=' ,F8.3,2X, ■PXA=' ,F8.3) 

82  CONTINUE 
C  PRINT  86,PXA,PYA 

C  86  FORMATC  '  ,  '  PXA=  •  ,  F8 . 3 , 2X,  '  PYA=  •,F8.3,/) 

C 

C  CALCULATE  THE  BASELINE  CORRECTION.  THERE  ARE  TWO  POSSIBLE 
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CASES:  (A)  POSITIVE  SLOPING  BASELINE,  OR  (B)  NEGATIVE 

SLOPING  BASELINE. 

CASE  (A)  -  POSITIVE  SLOPING  BASELINE  (BPl  <  BP2) 

SKIP  TO  CASE  (B)  IF  BASELINE  SLOPE  IS  NEGATIVE 
IF(BP1.GT.BP2)  GOTO  3000 

OTHERWISE  PROCEED  WITH  THE  CORRECTION 

DPl=(1309+1309)/2 

DP2= (1085+1080) /2 

YDIF=BP2-BPl 

XDIF=DP1-DP2 

XDIF2=DP1-PXA 

FAC1=XDIF2/XDIF 

FAC2=FAC1*YDIF 

FAC2=FAC2+BP1 

PRINT  87,  DP1,DP2, YDIF,XDIF,XDIF2,FAC1,FAC2 

87  FORMAT  (■  ‘  '  DPl=  *  ,  F8 . 3 , 2X,  * DP2=:  *  ,  F8 . 3  , 

2X, • YDIF=’ ,F9.4,2X, 

&'XDIF=‘ ,F9.4,2X, ■XDIF2= ‘ , F9 . 4 , 2X, / , •FAC1=' , 

&F9.4,2X, •FAC2=' ,F9.4,/) 

GOTO  3500 

CASE  (B)  -  NEGATIVE  SLOPING  BASELINE  (BPl  >  BP2) 

3000  CONTINUE 

DPl=(1309+1311)/2 
DP2= (1085+1080) /2 
YDIF=BP1-BP2 
XDIF=DP1-DP2 
XDIF2=PXA-DP2 
FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF 
FAC2=FAC2+BP2 
3500  CONTINUE 

APPLY  THE  CORRECTION  TO  THE  ABSORBANCE  VALUE  FROM  PPK 
CABS=PYA-FAC2 

PRINT  OUT  CORRECTED  ABSORBANCE  VALUE 
PRINT  88,  CABS 

88  FORMAT  ('  '  ,  ’CTFE::  •,F4.3) 

PRINT* , •  ' 

END  THE  SUBROUTINE 
RETURN 
END 

(NOTE:  this  subroutine  is  compiled  as  a  part  of  the  EXP4  batch  file.  See  APPENDIX  E 
for  the  listing] 
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APPENDIX  H 


SOFTWARE  LISTING  FOR  THE  FORTRAN  PROGRAM  USED  TO 
CALCULATE  THE  BASELINE-CORRECTED  ABSORBANCE 
VALUES  FOR  METHANE  VAPOR  AT  1305  CM-1 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c 

C  SUBROUTINE  CH4 (Y, FXAX.CABS) 

C 

C  CALCULATES  THE  BASELINE  BY  TAKING  THE  AVERAGE  OF  POINTS 
C  ON  ONE  SIDE  OF  THE  ANALYTICAL  PEAK  (3026  -  3023  CM-1) 

C  FOR  THE  BEGINNING  POINT.  THEN  TAKING  THE  AVERAGE  OF 
C  POINTS  ON  THE  OTHER  SIDE  OF  THE  ANALYTICAL  PEAK  (2996  - 
C  2992)  FOR  THE  ENDING  POINT.  A  LINE  CONNNECTING  THESE 
C  TWO  POINTS  FORMS  THE  BASELINE.  THE  ABSORBANCE  VALUE  AT 
C  THE  FIRST  ANALYTICAL  PEAK  (3018  CM-1)  IS  THEN  CORRECTED 
C  TO  THIS  MEASURED  BASELINE.  THE  SECOND  ANALYTICAL  PEAK 
C  AT  1305  CM-1  IS  CORRECTED  SIMILARLY.  (ONLY  THE  1305 
C  CM-1  PEAK  WAS  USED  IN  THIS  STUDY.) 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE  CH4(X,FXAX) 

REAL  X (22528) ,BP1,LBP,PXA,PYA,CABS,XDIF,XDIF2,YDIF, 
1FAC1,FAC2,EP,SP.IC 
INTEGER  COUNT 

DEFINE  THE  BASELINE  POINTS  FOR  THE  METHANE  BAND  AT  1305  CM-1 

IPT1=1310/FXAX 
IPT2=1308/FXAX 
IPT3=1280/FXAX 
IPT4=1279/FXAX 

PRINT  35,  IPT1,IPT2,IPT3,IPT4 

35  FORMAT ('  ' , ' IPTl= ' , I4 , 2X, ’ IPT2= • , I4 , 2X,  ' IPT3= ' , I4 , 2X, 
l‘IPT4=' ,14, /) 

AVERAGE  THE  POINTS  FOR  THE  FIRST  BASELINE  VALUE 

SP  =  0 
BP1=  0 
COUNT  =  0 
IFIRST=IPT1 
ILAST  =IPT2 

NOW  CALCULATE  FIRST  BASELINE  POINT 
PRINT  38,  IFIRST, ILAST 

38  FORMAT ('  ' , ' IFIRST= ' , 14 , 2X, ' ILAST= ' , I4 , / ) 

DO  72  I=ILAST, IFIRST 
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SP=SP+X(I) 

C0UNT=:C0UNT+1 

C 

C  CALCULATE  WAVENUMBER  VALUES  FOR  PRINTOUT 
C 

C  WN=(I-1) *FXAX 

C  PRINT  39, WN,X( I) ,SP, COUNT 

C  39  FORMATC  *  ,  ' WN=  ‘  , F8 . 3 , 2X, 'X  (I)  =  *  ,  F8 . 3 , 2X, 

C  &’SP=' ,F8.3,2X, 

C  l‘COUNT=' ,16) 

72  CONTINUE 

BP1=SP/C0UNT 

AVERAGE  THE  POINTS  FOR  THE  2ND  BASELINE  VALUE 

EP  =  0 
LBP  =  0 
COUNT  =  0 
IFIRST=IPT3 
ILAST  =IPT4 

NOW  CALCULATE  THE  SECOND  BASELINE  POINT 

DO  83  I=ILAST, IFIRST 
EP=EP+X ( I ) 

COUNT=COUNT+l 
83  CONTINUE 

LBP=EP/COUNT 

PRINT  OUT  THE  BASELINE  DETERMINING  POINTS 
PRINT  67,  BPl,LBP 

67  FORMATC  ■,'BPl=  '  ,  F9 . 4 , 2X,  '  LBP=  ■,F9.4,/) 

DETERMINE  THE  X,y  VALUES  (CM-1  AND  ABS . )  FOR  THE 
ANALYTICAL  PEAK  (1305  CM-1)  BY  FINDING  THE  MAXIMUM  Y  VALUE 
IN  THE  INTERVAL  FOR  THIS  PEAK  WHICH  IS;  1307  CM-1  TO  1304 
CM-1  (OR  IN  DATA  POINT  VALUES;  2710  TO  2704). 


PXA=0 

PYA=0 

DO  84  1=2704,2710 

IF  (X(I) .GT.X(I-l) )  THEN 
PYA=X ( I ) 

IC=I-1 

PXA=IC*FXAX 

ENDIF 

PRINT  399, I,X(I) ,PYA,IC,PXA 

399  FORMATC  '  ,  '  1=  '  ,  I4 , 2X,  •X= '  ,  F8 . 3 , 2X,  '  PYA=  ‘  ,  F8 . 3 , 2X, 

l'IC=' ,F8.3,2X, •PXA=’ ,F8.3) 

84  CONTINUE 

PRINT  89,PXA,PYA 

89  FORMATC  ’  ,  '  PXA=  '  ,  F8 . 3 , 2X,  ’  PYA=  •,F8.3,/) 


nnooooo  on  non  ooooo  non 


C  CALCULATE  THE  BASELINE  CORRECTION.  THERE  ARE  TWO  POSSIBLE 
C  CASES;  (A)  POSITIVE  SLOPING  BASELINE,  OR  (B)  NEGATIVE 
C  SLOPING  BASELINE. 

C 

C  CASE  (A)  -  POSITIVE  SLOPING  BASELINE  (BPl  <  LBP) 

C 

C  SKIP  TO  CASE  (B)  IF  BASELINE  SLOPE  IS  NEGATIVE 
C 

IF(BPl.GT.LBP)  GOTO  3100 

OTHERWISE  PROCEED  WITH  THE  CORRECTION 

DPl=(1310+1308)/2 
DP2= (1280+1279 ) /2 
yDIF=LBP-BPl 
XDIF=DPl-DP2 
XDIF2=DP1-PXA 
FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF 
FAC2=FAC2+BP1 

PRINT  187,  DP1,DP2, YDIF,XDIF,XDIF2,FAC1,FAC2 

187  FORMAT (•  • , '  DPl= ' , F8 . 3 , 2X, ' DP2= ■ , F8 . 3 , 2X, 

&‘YDIF=' ,F9.4,2X, 

&'XDIF=:'  ,F9.4,2X,  'XDIF2='  ,F9.4,2X,/,  •FAC1=*  , 

&F9.4,2X, ■FAC2=' ,  F9.4,/) 

GOTO  4500 

CASE  (B)  -  NEGATIVE  SLOPING  BASELINE  (BPl  >  LBP) 

3100  CONTINUE 

DP1= (1310+1308) /2 
DP2= (1280+1279) /2 
YDIF=BP1-LBP 
XDIF=DPl-DP2 
XDIF2=PXA-DP2 
FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF 
FAC2=FAC2+LBP 
4500  CONTINUE 

APPLY  THE  CORRECTION  TO  THE  ABSORBANCE  VALUE  FROM  PPK 
CABS=PYA-FAC2 

PRINT  OUT  CORRECTED  ABSORBANCE  VALUE 
PRINT  188,  CABS 

188  FORMAT (■  PEAK  ABS (CH4 , 1305  CM-1 ) =  ’ , F6 . 3 , / , 35  (  ' - ' ) , 

&!/) 

END  THE  SUBROUTINE 

RETURN 
END 

[NOTE:  this  program  is  compiled  and  linked  as  a  subroutine  of  EXP4.  See  Appendix  E  for 
the  applicable  batch  file  listing.] 
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APPENDIX  I 


SOFTWARE  LISTING  FOR  THE  FORTRAN  PROGRAM  USED  TO 
CALCULATE  THE  BASELINE-CORRECTED  ABSORBANCE 
VALUE  FOR  WATER  VAPOR  AT  1700  CM-1 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

C  SUBROUTINE  H20 ( Y, FXAX, CABS) 

C 

C  CALCULATES  THE  BASELINE  BY  TAKING  THE  AVERAGE  OF  POINTS 
C  ON  ONE  SIDE  OF  THE  ANALYTICAL  PEAK  (1703  -  1702  CM-1) 

C  FOR  THE  BEGINNING  POINT.  THEN  TAKING  THE  AVERAGE  OF 
C  POINTS  ON  THE  OTHER  SIDE  OF  THE  ANALYTICAL  PEAK  (1694  - 
C  1691  CM-1)  FOR  THE  ENDING  POINT.  A  LINE  CONNNECTING 

C  THESE  TWO  POINTS  FORMS  THE  BASELINE.  THE  ABSORBANCE 

C  VALUE  AT  THE  ANALYTICAL  PEAK  (1700  CM-1)  IS  THEN 
C  CORRECTED  TO  THIS  MEASURED  BASELINE. 

C 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

SUBROUTINE  H20 (Y, FXAX,CABS) 

REAL  Y (22528) , BPl,LBP, PXA, PYA,CABS ,XDIF, XDIF2 , YDIF, 
1FAC1,FAC2,EP,SP,IC 
INTEGER  COUNT 

DEFINE  THE  BASELINE  POINTS  FOR  THE  WATER  ANALYTICAL  PEAK 

IPT1=1703/FXAX 
IPT2=1702/FXAX 
IPT3=1694/FXAX 
IPT4=1691/FXAX 

PRINT  35,  IPT1,IPT2,IPT3,IPT4 

35  FORMATC  '  ,  '  IPTl=  *,14, 2X,  '  IPT2= ’,14 , 2X,  '  IPT3=  •  ,  I4 , 2X, 
l'IPT4=' ,14,/) 

AVERAGE  THE  POINTS  FOR  THE  1ST  BASELINE  VALUE 

SP  =  0 
BP1=:  0 
COUNT  =  0 
IFIRST=IPT1+1 
ILAST  =IPT2 

NOW  CALCULATE  FIRST  BASELINE  POINT 
PRINT  36,  IFIRST, ILAST 

36  FORMATC  '  ,  *  IFIRST=  *  ,  14 , 2X,  ’  ILAST=  *  ,  I4  ,  /  ) 

DO  70  I=ILAST, IFIRST 

SP=SP+Y(I+1) 

COUNT=COUNT+l 
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c 

C  CALCULATE  WAVENUMBER  VALUES  FOR  PRINTOUT 
C 

C  WN=(I)*FXAX 

C  PRINT  37,  WN, Y(I+1) ,SP, COUNT 

C  37  FORMATC  *  ,  '  WN=  '  ,  F8 . 3 , 2X,  ‘  Y  ( I )  =  '  ,  F8 . 3  , 

C  &2X, •SP=' ,F8.3, 

C  12X, •COUNT=* , 16) 

70  CONTINUE 

BPl=SP/COUNT 
C 

C  AVERAGE  THE  POINTS  FOR  THE  2ND  BASELINE  VALUE 
C 

EP  =  0 
LBP  =  0 
COUNT  =  0 
IFIRST=IPT3+1 
ILAST  =IPT4 
C 

C  NOW  CALCULATE  THE  SECOND  BASELINE  POINT 
C 

DO  80  I=ILAST, IFIRST 
EP=EP+Y(I+1) 

COUNT=COUNT+l 
C 

C  CALCULATE  WAVENUMBER  VALUES  FOR  PRINTOUT 
C 
C 
C 

C140 
C 

C137 
C 
C 

80 
C 

C  PRINT  OUT  THE  BASELINE  DETERMINING  POINTS 
C 

C  PRINT  85,  BP1,LBP 

C  85  FORMATC  '  ,  '  BPl=  '  ,  F9 . 4 , 2X,  ’  LBP=  •,F9.4,/) 

C 

C  DETERMINE  THE  X,Y  VALUES  (CM-1  AND  ABS .  )  FOR  THE 
C  ANALYTICAL  PEAK  (1700  CM-1)  BY  FINDING  THE  MAXIMUM  Y  VALUE 
C  IN  THE  INTERVAL  FOR  THIS  PEAK  WHICH  IS:  1702  CM-1  TO  1699 
C  CM-1  (OR  IN  DATA  POINT  VALUES:  3530  TO  3523). 

C 

PXA=0 

PYA=0 

DO  82  1=3523,3530 

IF  (Y(I) .GT.Y(I-l) )  THEN 
PYA=Y ( I ) 

IC=I-1 

PXA=IC*FXAX 


WN=(I) *FXAX 
PRINT  140 
FORMAT  (  '  '  ) 

PRINT  137,  WN,Y(I+1) ,SP, COUNT 
FORMATC  '  ,  •WN=‘  ,F8.3,2X,  '  Y  (I)  =  ’  ,  F8 . 3  , 
&2X, •SP=' ,F8.3, 

12X, ■COUNT=' , 16) 

CONTINUE 

LBP=EP/COUNT 
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ENDIF 

C  PRINT  299, I,Y(I) ,PYA,IC,PXA 

C299  FORMAT (•  * , *  1= • , I4, 2X, ‘ Y= • , F8 . 3 , 2X,  ‘ PYA= • , F8 . 3 , 2X, 

C  l’IC=' ,F8.3,2X, 'PXA=' ,F8.3) 

82  CONTINUE 

C  PRINT  86,PXA,PYA 

C  86  FORMAT (•  ‘ ' PXA=  ‘ , F8 . 3 , 2X, * PYA=  •,F8.3,/) 

C 

C  CALCULATE  THE  BASELINE  CORRECTION.  THERE  ARE  TWO  POSSIBLE 
C  CASES:  (A)  POSITIVE  SLOPING  BASELINE,  OR  (B)  NEGATIVE 
C  SLOPING  BASELINE. 

C 

C  CASE  (A)  -  POSITIVE  SLOPING  BASELINE  (BPl  <  LBP) 

C 

C  SKIP  TO  CASE  (B)  IF  BASELINE  SLOPE  IS  NEGATIVE 
C 

IF(BPl.GT.LBP)  GOTO  3000 

OTHERWISE  PROCEED  WITH  THE  CORRECTION 

DPl=(1703+1702)/2 
DP2= (1694+1691) /2 
YDIF=LBP-BP1 
XDIF=DP1-DP2 
XDIF2=DP1-PXA 
FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF 
FAC2=FAC2+BP1 

PRINT  87,  DPl,DP2,YDIF,XDIF,XDIF2,FACl,FAC2 
87  FORMATC  '  ,  'DPl=  '  ,  F8 . 3 , 2X,  '  DP2= '  ,  F8 . 3 , 

&2X, • YDIF=' ,F9.4,2X, 

&'XDIF=' ,F9.4,2X, 'XDIF2= • , F9 . 4, 2X, / , 'FAC1=' , 

&F9.4,2X, •FAC2=' ,F9.4,/) 

GOTO  3500 

CASE  (B)  -  NEGATIVE  SLOPING  BASELINE  (BPl  >  LBP) 

3000  CONTINUE 

DP1= (1703+1702) /2 
DP2= (1694+1691) /2 
YDIF=BP1-LBP 
XDIF=DP1-DP2 
XDIF2=PXA-DP2 
FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF 
FAC2=FAC2+LBP 

PRINT  57,  DP1,DP2,YDIF,XDIF,XDIF2,FAC1,FAC2 
57  FORMATC  ’  ,  'DP1=  ' , F8 . 3 , 2X,  •  DP2=  ‘  ,  F8 . 3 , 

&2X, •YDIF=* ,F9.4,2X, 

&’XDIF=' ,F9.4,2X, ‘XDIF2=* ,F9.4,2X, /, ’FAC1=’ , 

&F9.4,2X, •FAC2=’ ,  F9.4,/) 

3500  CONTINUE 

APPLY  THE  CORRECTION  TO  THE  ABSORBANCE  VALUE  FROM  PPK 
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c 

CABS=PYA-FAC2 

PRINT  OUT  CORRECTED  ABSORBANCE  VALUE 
PRINT  88,  CABS 

88  FORMAT (•  * , * CABS= * , F4 . 3 , / ) 

END  THE  SUBROUTINE 
RETURN 
END 

[NOTE;  this  program  was  compiled  and  linked  as  a  subroutine  of  EXP4..  See  Appendix 
E  for  the  applicable  batch  file  listing.] 


51 


APPENDIX  J 


SOFTWARE  LISTING  FOR  THE  FORTRAN  PROGRAM  USED  TO 
CALCULATE  THE  BASELINE-CORRECTED  ABSORBANCE 
VALUE  FOR  OZONE  AT  1055  CM-1 


ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c  c 

C  TAKES  THE  SET  BASELINE  POINTS  AND  CALCULATES  A  BASELINE-  C 
C  CORRECTED  ABSORBANCE  VALUE  FOR  THE  03  ANALYTICAL  PEAK  AT  C 
C  1055  CM-1  THEN  PRINTS  OUT  THE  VALUE.  C 
C  C 
C  THIS  VERSION  OF  THE  SUBROUTINE  HAS  BEEN  SPECIFICALLY  C 
C  DESIGNED  TO  MEASURE  OZONE  IN  THE  PRESENCE  OF  CTFE.  THE  C 
C  INTERFEREING  CTFE  ABSORBANCE  IS  FIRST  MEASURED  AND  THEN  C 
C  SCALED  AND  SUBTRACTED  OUT.  THIS  PROCESS  IS  DESCRIBED  C 
C  BELOW.  DETERMINE  THE  SCALE  FACTOR  FOR  SUBTRACTING  C 
C  OUT  CTFE  IN  THE  FOLLOWING  WAY;  C 
C  (1)  SELECT  A  REGION  OF  THE  CTFE  SPECTRUM  WHICH  DOES  NOT  C 
C  OVERLAP  ANY  OZONE  ABSORPTION  AND  FIND  THE  CTFE  C 
C  ABSORPTION  MAXIMUM  IN  THAT  REGION  BY  THE  USUAL  C 
C  BASELINE  CORRECTION  TECHNIQUE.  (THE  BAND  SELECTED  IS  C 
C  THE  ONE  WHICH  EXTENDS  FROM  1111  TO  1148  CM-1  WITH  C 
C  A  MAXIMUM  AT  1130.578  CM-1.)  C 
C  (2)  SUBTRACT  OUT  THE  CTFE  CONTRIBUTION  FROM  THE  OZONE  C 
C  SPECTRUM  IN  THE  1037  -  1074  CM-1  REGION  BY  SCALING  THE  C 
C  MEASURED  CTFE  ABSORBANCE  WITH  THE  ABSORBANCE  STORED  IN  C 
C  THE  DATA  STATEMENTS  THEN  SUBTRACTING  OUT  THE  CTFE  C 
C  CONTRIBUTION  TO  THE  ABSORBANCE  IN  THE  REGION  OF  INTEREST  C 
C  ON  A  POINT-BY- POINT  BASIS.  C 
C  C 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

SUBROUTINE  03 ( Y , FXAX , CABS ) 

REAL  Y(22528) , BPl, BP2 , PXA, PYA,CABS,XDIF, XDIF2 , YDIF, 
IFACI , FAC2 , EP , SP , REF ( 100 ) , SCALE, TEM, PXM , PYM, MAX 
INTEGER  COUNT 

DATA  REF/. 128, .130, .132, .133, .136, .140, .144,  .148, .151, 
1.156,  .162,  .167,  .171  , .175, .180, .184,  .187, .190, .192, .195, 

2.196. . 197. .197. .197. .196. .194. .192. .189. .187. .185. .180, 

3.176. . 172. .168. .163. .158. .152. .148. .145. .141. .137. .134, 

4.132. . 129. .127. .124. .123. . 122 . . 120 . . 118,  .116, .114, .111, 
5.108,  .105,  .101, .099, .096, .092, .088, .085, .083, .080, .078, 
6.076, .075, .074, .073, .073, .073, .073, .073, .073, .072, .073, 
7.073, .073, .074, .074/ 

C 

PXA=0 

PYA=0  Il=1124/FXAX  I2=1134/FXAX  MAX=0 
C  PRINT*, •  ■ 

C  PRINT*, ********  find  the  CTFE  PEAK  ABS .  AT  1130  CM-lC 

C  PRINT*, ■  ■ 
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DO  820  1=11,12 

IF  (ABS (Y(I) ) .GT.MAX)  THEN 
MAX=Y(I)  PYA=MAX 
IC=I-1  PXA=IC*FXAX 

C  PRINT  399,1, Y (I) ,PYA.PXA 

C399  FORMAT (•  ’ , ' 1= • , 14 . 2X. ' Y= * . F8 . 3 , 2X, \ 

C  &*PYA=' ,F8.3,2X, 

C  1 • PXA= • , F8 . 3 ) 

ENDIF 

820  CONTINUE 

C  PRINT  860 , PXA, PYA 

C860  FORMAT (■  ' , • PXA=  ' , F8 . 3 , 2X, * PYA=  •,F8.3,/) 

C  NOW,  MAKE  THE  BASELINE  CORRECTION  FOR  THIS  ABSORBANCE 

C  BY  DEFINING  A  BASELINE  FROM  ACTUAL  SPECTRAL  FEATURES. 

C 

C  DEFINE  THE  BASELINE  POINTS  FOR  THE  UPPER  END  OF  THE 
C  BAND  (1141-1146  CM-1) 

C 

IPT3=1141/FXAX  IPT4=1146/FXAX 
C 

C  AVERAGE  THE  POINTS  FOR  THE  HIGHER  CM-1  BASELINE  LOCATION 
C 

SP  =  0 
BP1=  0 
COUNT  =  0 
IFIRST=IPT3 
ILAST  =IPT4 
C  PRINT*,'  ' 

C  PRINT*,  *****  find  THE  CTFE  1130  CM-1  PEAK  BASELINE  PTS ' 

C  PRINT*,'  ' 

C 

C  NOW  CALCULATE  FIRST  BASELINE  POINT 
C 

IF1=IFIRST*FXAX  IF2=ILAST*FXAX 
C  PRINT  46,  IF1,IF2 

C  46  FORMATC  '  ,  '  FIRST=  '  ,  14 , 2X,  •LAST= '  ,  14 ,  / ) 

DO  71  I=IFIRST, ILAST 

SP=SP+Y(I)  COUNT=COUNT+l 
C 

C  CALCULATE  WAVENUMBER  VALUES  FOR  PRINTOUT 
C 

C  WN=(I-1)*FXAX 

C  PRINT  47,WN,Y(I) ,SP,COUNT 

C  47  FORMATC  '  ,  ' WN=  '  , F8 . 3 , 2X,  '  Y(I)  =  '  ,  F8 . 3 , 2X,  ' SP=  '  ,  F8 . 3 , 

C  12X, 'COUNT=' ,16) 

71  CONTINUE 

BPl=SP/COUNT 

C  PRINT*, 'BP1=  ',BP1 

C  PRINT* , '  ' 

C 

C  DEFINE  THE  BASELINE  POINTS  FOR  THE  LOWER  END  OF 
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THE  BAND  (1113-1117  CM-1) 

IPT5=1111/FXAX  IPT6=1117/FXAX 


AVERAGE  THE  POINTS  FOR  THE  LOWER  CM-1  BASELINE  LOCATION 

SP  =  0 
BP2=  0 
COUNT  =  0 
IFIRST=IPT5 
ILAST  =IPT6 

NOW  CALCULATE  LOWER  BASELINE  POINT 

IF2=IFIRST*FXAX  IL2=ILAST*FXAX 
PRINT  38,  IF2,IL2 

38  FORMAT {•  * , • FIRST= ' , I4 , 2X, • LAST= * , 14 , / ) 

DO  72  I=IFIRST, ILAST 

SP=SP+Y(I)  COUNT=COUNT+l 

CALCULATE  WAVENUMBER  VALUES  FOR  PRINTOUT 

WN=(I-1) *FXAX 

PRINT  39, WN, Y(I) ,SP, COUNT 

39  FORMATC  •  ,  'WN=' ,F8.3,2X, ’YtDs*  ,F8.3, 

&2X, 'SPs* ,F8.3, 

12X, 'COUNTS' ,16) 

72  CONTINUE 

BP2=SP/C0UNT 
PRINT*, 'BP2=  ',BP2 
PRINT* , '  ' 

CALCULATE  THE  BASELINE  CORRECTION.  THERE  ARE  TWO  POSSIBLE 
CASES;  (A)  POSITIVE  SLOPING  BASELINE,  OR  (B)  NEGATIVE 
SLOPING  BASELINE. 

CASE  (A)  -  POSITIVE  SLOPING  BASELINE  (BPl  <  BP2) 

SKIP  TO  CASE  (B)  IF  BASELINE  SLOPE  IS  NEGATIVE 
IF(BP1.GT.BP2)  GOTO  3010  C  C  OTHERWISE  PROCEED  WITH  THE 
CORRECTION 

PRINT*,'  ' 

PRINT*,'****  CASE  (A)  BASELINE  CORRECTION  ***' 
PRINT*,'  (CTFE  REFERENCE  PEAK  AT  1130  CM-1' 

PRINT*, '  ' 

DP1= (1141+1146) /2  DP2= (1111+1117) /2  YDIF=BP2-BPl 
XDIF=DP1-DP2  XDIF2=DP1-1131  FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF  FAC2=FAC2+BPl 
PRINT  88,  DP1,DP2, YDIF,XDIF,XDIF2,FAC1,FAC2 
88  FORMATC  '  ,  ' DPl=  '  ,  F8 . 3 , 2X, 'DP2=  '  ,  F8 . 3 , 2X,  '  YDIF=  '  , 

&F9.4,2X,  'XDIF=' ,F9.4,2X, 'XDIF2=' ,F9.4,2X, /, 'FAC1=' , 
&F9.4,2X, 'FAC2=' ,F9.4,/) 
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no  o  o  o  o  o  o  n  o  o  o  o  o  o  o  o  o  o  o  o  n  o  o  o  o  o  w  o  o  o  o  o  o  o  o  oj  o  o  o 


GOTO  3510 


CASE  (B)  -  NEGATIVE  SLOPING  BASELINE  (BPl  >  BP2) 

010  CONTINUE 
PRINT*, •  • 

PRINT*,'****  CASE  (B)  BASELINE  CORRECTION  ****• 
PRINT*, •  (CTFE  REFERENCE  PEAK  AT  1130  CM-1) ' 

PRINT*,'  ' 

DP1= (1141+1146) /2  DP2= (1111+1117) /2  yDIF=BPl-BP2 
XDIF=DP1-DP2  XDIF2=1131-DP2  FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF  FAC2=FAC2+BP2 
PRINT  58,  DP1,DP2, YDIF,XDIF,XDIF2,FAC1,FAC2 
58  FORMATC  ‘  ‘ DP1= '  , F8 . 3 , 2X,  ' DP2= '  ,  F8 . 3 , 2X,  '  YDIF=  '  , 
&F9.4,2X, 'XDIF=’ ,F9.4,2X, •XDIF2='  ,F9. 4, 

&2X./, ■FAC1=' ,F9.4,2X, •FAC2=' ,F9.4, /) 

510  CONTINUE 

THE  BASELINE-CORRECTED  ABSORBANCE  VALUE  AT  1130.578  CM-1 
CAN  NOW  BE  USED  TO  SUBTRACT  OUT  THE  CTFE  FROM  THE 
OZONE  SPECTRUM 


THE  SCALE  FACTOR  IS  DETERMINED  BY  THE  RATIO  OF  THE 
BASELINE -CORRECTED  CTFE  ABSORBANCE  VALUE  AT  1130.578  CM-1 
IN  THE  REFERENCE  SPECTRUM  (MEASURED  TO  BE  0.270)  AND 
THE  BASELINE-CORRECTED  CTFE  ABSORBANCE  VALUE  AT  COMPUTED 
MAXIMUM  ABSORBANCE  IN  THE  SAMPLE  SPECTRUM 

AREF  =  PYA-FAC2 
PRINT  89,  AREF 

89  FORMAT(35 (•-'),/,  '  ',' PEAK  ABS (1130  CM-1) =  ',F6.3,/, 
135('-'),/) 

SCALE  =  ABS (AREF) / (0.270) 

PRINT  * , ' AREF=  ’ , AREF , ' SCALE= ' , SCALE 

SUBTRACT  OUT  THE  CTFE  CONTRIBUTION  TO  THE  ABSORBANCE 
FROM  1037  TO  1074  CM-1  SO  THAT  OZONE  ABSORBANCE  CAN  BE 
COMPUTED  ACCURATELY 
PRINT*, •  • 

PRINT*, '***  SUBTRACT  OUT  ANY  CTFE  IN  THE  03  BAND' 

PRINT*, ■  • 

IPT5  =  1037/FXAX 
IPT6  =  1074/FXAX 
DO  40  1=  IPT5,IPT6 
J=I+1-IPT5 
TEM=Y(I) 

REF ( J) =SCALE*REF ( J) 

Y(I)=Y(I)-REF(J) 

L=I*FXAX 

M=J*FXAX 

PRINT  42,L,TEM,Y(I) ,REF(J)  ,M 

42  FORMATC  '  ,  '  1= '  .  16 , 2X,  '  Y ( I)  =  *  , F6 . 3 , 2X,  '  Y ( I )  CORR.  =  '  , 
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1F6.3,2X, ‘REF(J)=‘ ,F6.3,2X. •J=’ ,14) 

0  CONTINUE 

DEFINE  THE  BASELINE  POINTS  FOR  THE  HIGHER  CM-1  OZONE 
ANALYTICAL  PEAK 

PRINT*, •  ’ 

PRINT*, *****  NOW,  PRINT  OUT  THE  BASELINE  CORRECTION* 
PRINT*, •  VALUES  FOR  THE  03  BAND' 

PRINT*,'  ■ 

IPT1=1070/FXAX  IPT2=1074/FXAX 
PRINT  45,  IPT1,IPT2 

45  FORMAT ('  ' , ' IPTl= ' , 14 , 2X, ' IPT2= '  ,  I4 ,  / ) 

AVERAGE  THE  POINTS  FOR  THE  HIGHER  CM-1  BASELINE  LOCATION 

SP  =  0 
BP1=  0 
COUNT  =  0 
IFIRST=IPT1 
ILAST  =IPT2 

NOW  CALCULATE  FIRST  BASELINE  POINT 

IF1=IFIRST*FXAX  ILl=ILAST*FXAX 
PRINT  36,  IFl.ILl 

36  FORMATC  ' ,  '  IFIRST=  '  ,  14 , 2X,  '  ILAST= ' ,  14, /) 

DO  70  I=IFIRST, ILAST 

SP=SP+Y(1)  C0UNT=C0UNT+1 

CALCULATE  WAVENUMBER  VALUES  FOR  PRINTOUT 

WN=(I-1)*FXAX 

PRINT  37,WN, Y(I) ,SP,COUNT 

37  FORMATC  '  ,  ' WN=  '  ,  F8 . 3 , 2X,  '  Y ( I)  =  '  ,  F8 . 3 , 2X, 

&'SP=' ,F8.3, 

12X, 'COUNT=' ,16) 

70  CONTINUE 

BP1=SP/C0UNT 
PRINT*, 'BP1=  •,BP1 

THE  LOWER  CM-1  BASELINE  POINT  IS  DEFINED  AS  THE  VALLEY 
MINIMUM 

IN  THE  RANGE  1042-1045  CM-1  (OR  IN  DATA  POINT  VALUES: 
2161-2167).  THIS  SECTION  FINDS  THE  MINIMUM  VALUE  FOR  THE 
ABSORBANCE  AND  RECORDS  THE  CORRESPONDING  CM-1  VALUE. 

PRINT*, ■  ■ 

PRINT*, *****  find  the  03  PEAK  MIN.  IN  THE  RANGE  OF* 
PRINT*,'  1042-1045  CM-1 • 

PRINT*,  *  * 

PXM=  0 
PYM  =  0 
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DO  75  K=2161,2167 

IF  (Y(K) .LT.Y(K-l) )  THEN 

PYM=Y(K)  IE=K-1  PXM=IE*FXAX 
C  PRINT  199,K,Y(K) ,PYM,PXM 

C199  FORMAT (■  • . ‘ K= ' , I4 , 2X, * Y= ' , F8 . 3 , 2X, 

C  &'PYM=' ,F8.3.2X, •PXM=‘ ,F8.3) 

ENDIF 

75  CONTINUE 

BP2=PYM 

C  PRINT*, •BP2=  ■,BP2 

C 

C  DETERMINE  THE  X.Y  VALUES  (CM-1  AND  ABS.)  FOR  THE  OZONE 
C  ANALYTICAL  PEAK  (1055  CM-1)  BY  FINDING  THE  MAXIMUM  Y  VALUE 
C  IN  THE  INTERVAL  FOR  THIS  PEAK  WHICH  IS:  1052  CM-1  TO  1058 
C  CM-1  (OR  IN  DATA  POINT  VALUES:  2182  TO  2194)  . 

C 

C 

C  PRINT*, ‘  • 

C  PRINT*,'****  FIND  THE  03  PEAK  ABS.  AT  1055  CM-1' 

C  PRINT*, '  ' 

PXA=0 

PYA=0 

MAX=0 

DO  82  1=2182,2194 

IF  (ABS(Y(I) ) .GT.MAX)  THEN 
MAX=Y(I)  PYA=MAX 
IC=I-1  PXA=IC*FXAX 

C  PRINT  299,1, Y (I) ,PYA,PXA 

C299  FORMATC  '  ,  '  1=  '  ,  I4 , 2X,  '  Y=  '  , F8 . 3 , 2X,  '  PYA=  '  , 

C  &F8.3,2X, 'PXA=' ,F8.3) 

ENDIF 

82  CONTINUE 

C  PRINT  86,PXA,PYA 

C  86  FORMATC  '  ,  '  PXA=  '  ,  F8 . 3 , 2X,  '  PYA=  ■,F8.3,/) 

C 

C  CALCULATE  THE  BASELINE  CORRECTION.  THERE  ARE  TWO  POSSIBLE 
C  CASES:  (A)  POSITIVE  SLOPING  BASELINE,  OR  (B)  NEGATIVE 
C  SLOPING  BASELINE. 

C 

C  CASE  (A)  -  POSITIVE  SLOPING  BASELINE  (BPl  <  BP2 ) 

C 

C  SKIP  TO  CASE  (B)  IF  BASELINE  SLOPE  IS  NEGATIVE 
C 

IF{BP1.GT.BP2)  GOTO  3000 
OTHERWISE  PROCEED  WITH  THE  CORRECTION 
PRINT*,'  ' 

PRINT*,'***  CASE  (A)  BASELINE  CORRECTION  FOR  03' 
PRINT* , '  ' 

DP1= (1070+1074) /2  DP2= (1042+1045) /2  YDIF=BP2-BPl 
XDIF=DP1-DP2  XDIF2=DP1-PXA  FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF  FAC2=FAC2+BPl 
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PRINT  87,  DP1,DP2, YDIF,XDIF,XDIF2,FAC1,FAC2 
87  FORMAT ('  ‘ ' DPl= • , F8 . 3 , 2X, * DP2= ‘ , F8 . 3 , 2X, ‘ yDIF= • , 
&F9.4,2X, •XDIF=' ,F9.4,2X, •XDIF2=* ,F9.4, 

&2X,/, •FAC1=' ,F9.4.2X. •FAC2= * , F9 . 4 , / ) 

GOTO  3500 

CASE  (B)  -  NEGATIVE  SLOPING  BASELINE  (BPl  >  BP2 ) 

000  CONTINUE 


PRINT*,*  * 

PRINT*,'***  CASE  (B)  BASELINE  CORRECTION  FOR  03* 

PRINT*, '  • 

DP1= (1070+1074) /2  DP2= (1042+1045) /2  YDIF=BPl-BP2 
XDIF=DP1-DP2  XDIF2=PXA-DP2  FAC1=XDIF2/XDIF 
FAC2=FAC1*YDIF  FAC2=FAC2+BP2 
PRINT  57,  DP1,DP2, YDIF,XDIF,XDIF2,FAC1,FAC2 
57  FORMAT (•  * , * DPl= ' , F8 . 3 , 2X, * DP2= * , F8 . 3 , 2X, ' YDIF= * , 
&F9.4,2X, ■XDIF=' ,F9.4,2X, •XDIF2=' ,F9.4, 

&2X,/, •FAC1=' ,F9.4,2X, * FAC2= * , F9 . 4 , / ) 

3500  CONTINUE 

APPLY  THE  CORRECTION  TO  THE  ABSORBANCE  VALUE  FROM  PPK 
CABS=PYA-FAC2 

PRINT  OUT  CORRECTED  ABSORBANCE  VALUE 
PRINT  90,  CABS 

90  FORMAT(35( '-'),/, *  ',* PEAK  ABS (03 , 1055  CM-1 ) =  ',F6.3,/, 
135(*-*),/) 


END  THE  SUBROUTINE 
RETURN  END 
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APPENDIX  K 


LISTING  OF  MS-EXCEL  MACROS  USED  FOR  DATA  FORMATING 
(see  APPENDIX  A  for  a  description  of  how  the  ASCII  data  file  is  acquired 

and  sent  to  MS-EXCEL) 


OPERATION:  Once  the  appropriate  ASCII  data  file  is  transferred  to  MS-EXCEL,  the  data  is 
parced  manually  and  the  top  4  lines  of  the  worksheet  are  left  blank.  Then  the  two  macros  listed 
here  are  run.  Finally,  the  data  columns  are  labeled  according  to  the  chemical  species  which  are 
present  in  the  experiments  being  conducted,  lire  macros  are  stored  under  the  file  name 
NICOLET  PROCESSING  MACROS. 


ILA 


Recordl 

=SELECT("C1:C15") 

=C0PV()  ^ 

=PAS'rE() 

=SELECT("C1:C4") 

=INSeRT(l) 

=SfeLCCt("Rl") 

=SELEciT(Vi:ti2") 

=AL1<jNMENT(3) 

^LEcrrcn 

=EORMAT.NUMBER("d-mmm-yy") 

=SELECT("C2") 

=FORMAT.broMBER(''0.0000'“) 

=SELECT("C4") 

=E6RMAT.broMBEk("0.00<y0'^) 

=SELEcT("C5") 

=FORMAT.NUMbER{''(5EbtERAL“) 

=SELfiCT("C^‘') 

”FORMAT.NtJM6£R(''0.0") 

=SELECTrR2C4"i 

=EORKltJLA^"=ftCi[21+kC[31/60+RC[4]/3600’') 

=sele(::!t("R2C3") 

=FORMt]LAr=lF(Rt- 
1 1C[  1 1>RC[  1 1,IF(RC[  1 1>0.000 1 , 1 ,0),0)" ) 

=SELECiT(’^R2C2") 

=FoRMuLAr=i4*sUM(R2C5:kCtTr^+k(:ir51- 

R2C4'*) 

=SELECT("R2Cr) 

=FORMl)LA("=RC[4|+SUM?R2C3:RCr21)‘'5 

59 


“  rti'vv 


B.t:l«:f»irtf.lF:i _ 

B2£il^iSSl£!rMlli£SKl%tWl 


lEiW#;T.Tireiii:inai 

-!.iaii:c.iifgte]r.iH 


&1  i  :CN  ai* 

B>1^y:C«iir4iWIWfci 

EcmWgiNifMTIPr^i 


E»]aaHIKTIiW;l.-«iiiar*l!l;^ 


•2slSs3IUSI2!S£l 


lEwiiWitiKTiiii.tiiriiiaw 


IEt»l;WHIKTIilMBIJW\ 

mmmwm 


6 


